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ABSTRACT 


The  dynamic  response  of  the  oceanic  mixed-layer  to  the  thennodynamic  forcing 
at  the  sea  surface  is  analyzed  in  order  to  describe  the  pattern  of  the  oceanic  meridional 
overturning.  The  technique  proposed  in  this  study  is  based  on  residual-mean  theory, 
which  takes  into  account  the  transport  of  buoyancy  and  tracers  by  transient  eddies.  From 
the  observed  air-sea  density  flux  and  mixed-layer  density  distributions,  we  estimate  the 
two  components  of  the  Meridional  Overturning  Circulation  (MOC)  corresponding  to  the 
adiabatic  (along-isopycnal)  advection  and  the  diabatic  (cross-isopycnal)  flux. 
Calculations  are  performed  for  the  global  ocean  and,  additionally,  for  each  oceanic  basin. 
The  derived  method  extends  the  Walin  (1982)  water  mass  transformation  theory,  and 
pennits,  for  the  first  time,  assessment  of  the  strength  of  the  MOC  adiabatic  component 
from  the  sea  surface  data. 

This  study  offers  a  statistical  description  of  the  atmospheric  and  oceanic 
databases  and  gives  some  suggestions  for  the  choice  of  specific  datasets.  In  particular,  the 
two  most  reliable  atmospheric  climatology  databases  (ECMWF  and  NCEP/NCAR  re¬ 
analyses)  are  compared,  and  the  impact  of  their  inaccuracies  on  the  MOC  calculations  is 
evaluated. 

The  results  presented  in  this  paper  are  consistent  with  the  pattern  of  the 
thermohaline  circulation  estimated  from  in-situ  measurements  and  models.  They  also 
support  the  previous  estimates  of  diapycnal  volume  flux  by  the  Walin-type  calculations. 
Furthermore,  they  suggest  that  the  global  thennohaline  circulation  is  dominated  by  the 
adiabatic  advection  in  the  ocean  interior. 

One  of  the  goals  of  the  data  analysis  and  numerical  modeling  effort  is  to  explain 
the  role  of  the  mesoscale  variability  in  the  dynamic  coupling  of  the  ocean  and  the 
atmosphere.  Properties  of  the  ocean  thermal  fronts  and  eddies  detennine  the  Undersea 
Warfare  (USW)  tactics  in  the  areas  of  high  mesoscale  activity,  and  therefore  our  efforts 
to  predict  their  distribution  and  strength  are  directly  related  to  the  Navy  research  interests. 
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I.  INTRODUCTION 


The  ocean  plays  a  major  role  in  the  climate  system  since  it  is  the  main  absorber  of 
solar  energy  and  therefore  constitutes  the  direct  source  of  most  of  the  energy  in  the 
atmosphere.  Thus,  a  complete  understanding  of  the  complex  interactions  and  feedbacks 
between  the  ocean  and  the  atmosphere  is  essential  for  forecasting  any  climate  variation. 

The  interconnections  between  the  heat  and  moisture  fluxes  at  the  sea  surface  and 
the  oceanic  and  atmospheric  circulations  detennine  the  equilibrium  state  of  the  climate 
system.  Of  particular  importance  are  the  meridional  transports  of  energy  by  winds  and 
currents  (wind-  and  density-driven).  As  a  consequence,  fluctuations  in  the  density-driven 
oceanic  Meridional  Overturning  Circulation  (MOC),  also  referred  to  as  thermohaline 
circulation,  may  lead  to  major  climate  variations.  The  MOC  has  recently  become  a 
subject  of  great  interest  in  the  context  of  the  currently  observed  global  warming  in  polar 
regions. 

One  of  the  biggest  uncertainties  in  the  theory  of  the  MOC  is  associated  with  the 
relative  importance  of  the  interior  mixing  and  the  adiabatic  advection  for  its 
establishment  and  maintenance.  The  classical  theories  use  purely  diffusive  models  to 
depict  the  thermohaline  circulation:  Stommel  (1961)  emphasized  the  importance  of  the 
thennodynamic  forcing  on  the  individual  density  components,  and  Munk  (1966) 
described  the  vertical  mixing  as  the  means  for  communicating  the  buoyancy  anomaly 
from  the  surface  into  the  ocean  interior.  Even  more  recent  comprehensive  discussions 
(Welander  1986,  and  Whitehead  1995)  upheld  such  a  description  of  the  MOC. 

However,  several  new  theoretical  and  observational  findings  cast  some  doubt  on 
the  ability  of  these  purely  diffusive  models  to  explain  the  magnitude  and  pattern  of  the 
meridional  overturning.  Indeed,  although  numerical  simulations  of  the  Atlantic  exhibit  a 
dependence  of  the  MOC  on  diapycnal  diffusivity  (Bryan  1987,  among  others),  most 
diffusive  models  have  to  rely  on  very  high  values  of  vertical  mixing  (5-10x10 ~5m2s~l ) 
to  reproduce  oceanic  overturning  of  realistic  strength.  Yet  these  high  values  are  not 
supported  by  the  microstructure  and  tracer  dispersion  measurements  which  instead 
suggest  the  thermocline  diffusivity  on  the  order  of  1x10  5m2s  1  (Ledwcll  et  al.  1993; 
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Toole  et  al.  1994).  This  has  led  to  an  alternative  adiabatic  view  of  the  meridional 
overturning,  which  emphasizes  the  importance  of  the  advection  along  isopycnal  surfaces, 
induced  by  wind  forcing  and  eddy  transfer  (Toggvieler  and  Samuels  1998;  Gnanadesikan 
1999;  Marshall  and  Radko  2003;  Radko  2006).  The  relative  importance  of  the  interior 
mixing  and  adiabatic  advection  for  the  MOC  establishment  and  maintenance  has  become 
the  subject  of  a  vigorous  ongoing  debate. 

The  classical  representation  of  the  MOC,  shown  in  Fig.  1  and  Fig.  2,  focuses  on 
the  abyssal  circulation  driven  by  a  few  deepwater  sources  where  deep  convection  occurs, 
and  a  large-scale  upwelling  that  compensates  for  the  localized  downward  flux  (Stommel 
1961).  The  modem  theory  of  meridional  overturning  (Webb  and  Suginohara  2001; 
Bocaletti  et  al.  2005;  among  others)  identifies  and  describes  at  least  two  distinct  dynamic 
components  of  circulation,  which  are  the  shallow  overturning  cells  in  the  main 
thennocline  and  the  deep  circulation  in  the  abyssal  ocean.  In  the  abyssal  regions,  shielded 
from  the  direct  influence  of  wind,  small-scale  mixing  processes  are  necessary  to  re¬ 
supply  the  potential  energy  removed  in  the  interior  by  the  overturning  and  eddy¬ 
generating  process  (Wunsh  and  Ferrari  2004).  On  the  other  hand,  the  dynamics  of  the 
strongly  stratified  thennocline  and  the  associated  shallow  overturning  are  controlled  by 
the  ventilation  of  water  masses  along  the  isopycnals  that  outcrop  at  the  sea-surface 
(Luyten  et  al.,  1983);  therefore,  the  role  of  diapycnal  mixing  in  the  maintenance  of  the 
upper  cell  is  questionable  (Marshall  et  al.  2002).  Indeed,  unlike  the  abyssal  regions,  the 
thennocline  is  directly  influenced  by  the  surface  winds,  which  provide  the  energy  for  the 
meridional  overturning  and  make  it  possible  to  close  the  circulation  even  in  the  absence 
of  vertical  diffusion  below  the  mixed-layer. 

Traditional  views  on  the  oceanic  overturning  tended  to  emphasize  the  dynamics 

of  the  deep  overturning  circulation  (Roemmich  and  Wunsh  1985;  Munk  and  Wunsh 

1998)  and  its  major  role  in  the  meridional  heat  transport.  However,  the  relative 

significance  of  the  upper  and  deep  circulation  cells  should  be  re-evaluated  in  the  context 

of  the  recent  suggestion  (Boccaletti  et  al.  2005)  that  the  meridional  transport  of  heat  is 

controlled  by  the  processes  operating  in  the  upper  ocean.  In  an  attempt  to  meaningfully 

quantify  partitioning  of  the  oceanic  heat  flux  between  the  deep  and  shallow  overturning 

cells,  these  authors  introduced  the  “heatfunction”  which  identifies  the  components  of 
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circulation  that  effectively  contribute  to  the  total  meridional  heat  flux.  As  shown  in  Fig. 
3,  diagnostics  of  a  numerical  model  (Boccaletti  et  al.  2005)  have  pointed  out  that  the 
heatfunction  maxima  are  associated  with  a  surface  intensified  flow  largely  limited  to  the 
central  thermocline.  This  suggests  that  the  meridional  heat  transport  is  dominated  by  the 
contribution  from  the  upper  branch  of  the  MOC,  the  dynamics  of  which  are  still  poorly 
understood  (Samelson  2004). 


heat  loss  ATMOSPHERE  heat  input 

from  ocean  to  ocean 


1  l  i 


Figure  1.  Schematic  of  the  MOC  (from  http://www.Iiv.ac.uk/physocean/research.htm). 


Although  the  global  characteristics  of  the  oceanic  overturning  can  be  estimated 
from  the  air-sea  fluxes  (Trenberth  and  Caron  2001),  the  vertical  pattern  of  the  MOC  is 
difficult  to  reconstruct  from  observations  or  models.  Sparse  oceanographic  measurements 
allow  only  a  crude  partitioning  of  the  MOC  in  terms  of  the  volume  fluxes  in  the  shallow, 
intennediate,  and  deep  waters  (Talley  2003).  All  numerical  modeling  studies  rely  on  the 
assumed  parameterizations  of  diapycnal  fluxes,  and  the  magnitude  and  functional 
dependencies  of  these  fluxes  are  still  a  greatest  source  of  uncertainty.  This  difficulty,  for 
instance,  precludes  a  reliable  quantitative  estimate  of  the  relative  importance  of  the 
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adiabatic  advection  and  cross-isopycnal  diffusion  for  maintenance  of  the  MOC.  A 
successful  approach  in  analyzing  the  meridional  overturning  was  proposed  by  Walin 
(1982),  who  deduced  the  integrated  rates  of  water  mass  transformation  at  isopycnal 
surfaces  from  the  air-sea  buoyancy  fluxes  and  surface  density  data.  Walin’s  technique 
was  extended  and  applied  to  the  Atlantic  (Speer  and  Tziperman  1992;  Donners  et  al. 
2005)  and  the  global  ocean  (Speer  et  al.  1995).  The  water  mass  transformation  at  a  given 
buoyancy  surface  is  directly  linked  to  the  cross-isopycnal  mass  flux  (Tziperman,  1986), 
and  this  provides  an  estimate  of  the  strength  of  the  diapycnal  (and  therefore  inherently 
diabatic)  component  of  the  oceanic  circulation. 


A 


samw  Subantarctic  Mode  Water 

AAIW  Antarctic  Intermediate  Water 

RSOW  Red  Sea  Overflow  Water 

AABW  Antarctic  Bottom  Water 

NPDW  North  Pacific  Deep  Water 

AAC  Antarctic  Circumpolar  Current 

CDW  Circumpolar  Deep  Water 

NADW  North  Atlantic  Deep  Water 

UPPER  IW  Upper  Intermediate  Water 

10 DW  Indian  Ocean  Deep  Water 


Figure  2.  Classical  representation  of  the  global  thermohaline  circulation  and  its  associated 
oceanic  conveyor  belt  (from  Aguado  and  Burt,  Understanding  Weather  and 
Climate). 
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Assessing  the  adiabatic  pole-to-pole  component  of  the  MOC  is  more  challenging 
since  the  movement  of  fluid  along  the  isopycnal  surfaces  does  not  significantly  involve 
the  interior  water  mass  transformation.  Another  complication  results  from  the  presence 
of  mesoscale  eddies  whose  contribution  to  the  overall  overturning  has  to  be  taken  into 
account.  It  is  now  widely  accepted  that  the  eddy  transfer  of  buoyancy  and  tracers  in  the 
upper  ocean  is  often  comparable  in  magnitude  to  the  advection  by  the  time-mean  flow 
(e.g.  Radko  and  Marshall  2004;  Henning  and  Vallis  2004).  Since  eddies  generally  tend 
to  counteract  the  Eulerian  transport,  the  residual  flow  (sum  of  the  mean  and  eddy-induced 
circulations)  is  considerably  different  from  the  mean.  Nevertheless,  attempts  have  been 
made  to  diagnose  the  residual  overturning  circulation  from  the  air-sea  fluxes  and  surface 
buoyancy  distribution  (Speer  et  al.  2000;  Karsten  and  Marshall  2002;  Marshall  and 
Radko  2003).  These  studies  focused  on  the  southern  ocean,  where  predominantly  zonal 
flow  of  the  Antarctic  Circumpolar  Current  justified  the  use  of  a  simplified  two- 
dimensional  framework.  A  clear  way  of  diagnosing  the  vertical  structure  of  the  global 
pole-to-pole  MOC  from  the  air-sea  fluxes  for  complicated  three-dimensional  geometries 
of  the  isopycnals  is  still  lacking. 

The  purpose  of  this  study  is  to  combine  Walin’s  technique  with  the  ideas  of  the 
residual  circulation  theory  to  quantify  both  the  diabatic  and  the  adiabatic  components  of 
the  meridional  overturning.  More  precisely,  the  along-  and  cross-isopycnal  residual 
fluxes  in  the  whole  ocean  and  in  each  basin  (Atlantic,  Pacific,  and  Indian)  will  be 
estimated  from  the  air-sea  fluxes  and  mixed-layer  buoyancy  distribution.  The  technique 
presented  here  offers  for  the  first  time  a  quantitative  assessment  of  the  adiabatic 
component  of  the  MOC  which  is  based  entirely  on  the  sea  surface  data. 

The  coordinate  system  introduced  in  this  study  is  somewhat  similar  to  the 
isentropic  framework  often  used  by  meteorologists.  Here,  however,  the  isentrope  field 
cannot  represent  the  density  distribution,  since  the  seawater  potential  density  also 
depends  significantly  on  salinity.  In  addition,  contrary  to  the  isentropic  circulation,  the 
thermodynamic  forcing  in  this  MOC  model  is  restricted  to  the  sea-surface  and  the 
associated  oceanic  boundary  layer.  The  circulation  in  the  ocean  interior  is  assumed  to  be 
a  linear  response  to  mixed-layer  dynamics  induced  by  the  air-sea  fluxes. 
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The  formulation  of  the  diagnostic  technique  is  presented  in  section  II.  Section  III 
compares  the  different  databases  and  suggests  the  choice  of  specific  datasets.  In  section 
IV,  the  diagnostic  model  is  applied  to  the  global  ocean  and  separately  to  the  Atlantic, 
Pacific,  and  combined  Indian-Pacific  basins.  Finally,  the  main  results  and  the  MOC 
characteristics  are  discussed  in  the  conclusion  (section  V). 
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Figure  3.  (a)  Stream  function  related  to  the  zonally-averaged  flow  inferred  from  the  M.I.T. 

global  circulation  model.  The  distribution  of  meridional  volume  transport,  given 
in  Sv  ( 1 06  1 ),  is  defined  by  the  stream  function  since  seawater  is  nearly 

incompressible.  Red  shading  corresponds  to  clockwise  flow  and  blue  to  counter¬ 
clockwise  flow.  Contour  interval  is  5  Sv.  (b)  Zonal  averages  of  the  total  heat 
transport  of  the  global  circulation  model,  (c)  Meridional  heat  function  in  PW  (1015 
W)  inferred  from  the  global  circulation  model.  Red  shading  corresponds  to 
northward  heat  transport  and  blue  to  southward  heat  transport.  Contour  interval  is 
0.17  PW  (from  Boccaletti  et  al.  2005). 
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II.  FORMULATION 


The  two  surface  properties  used  to  reconstruct  the  large  scale  meridional 
overturning  in  the  ocean  are  the  mixed-layer  density  anomaly  ( <jm )  and  the  air-sea 

density  flux  ( B  ).  The  dimensionless  mixed-layer  density  anomaly  ( crm  )  is  related  to  the 

mixed-layer  density  pm  (in  kg  m  3 )  as  follows: 


Pm 


=  1000  +  o-,„ 


And  the  sea  surface  density  flux  ( B  )  is  conventionally  defined  as: 


D  ccH  a  (E-P)S 
B  = - +  M“ - 


C„ 


1-5 


(1) 


where  (a,/?)  are  the  expansion/contraction  coefficients  of  seawater,  p0  is  the  standard 
seawater  density,  H  is  the  heat  flux  from  the  atmosphere  into  the  ocean,  E  and  P  are  the 
evaporation  and  precipitation  rates,  S  is  the  salinity,  and  Cp  is  the  specific  heat  capacity 

of  water.  The  technique  utilizes  the  observational  data  based  on  meteorological  re¬ 
analysis  (density  flux  at  the  sea  surface)  and  Levitus  database  (seawater  density)  to 
reconstruct  the  pattern  of  the  meridional  overturning. 

For  the  computations  in  the  Atlantic,  Pacific,  and  combined  Indian-Pacific 
oceans,  we  consider  the  density  range  where  all  isopycnals  are  zonally  bounded  by  land, 
as  sketched  in  Fig.  4.  Confining  the  analysis  to  these  water  masses  greatly  simplifies  the 
interpretation  of  the  diagnostic  results  and  accurately  describes  the  upper  ocean.  For  the 
global  ocean,  computations  are  extended  to  the  maximum  value  of  mixed-layer  density 
anomaly. 


A.  RESIDUAL  FLUX  AT  THE  BOTTOM  OF  THE  MIXED-LAYER 

It  has  been  shown  (Andrews  and  McIntyre  1976)  that  the  distribution  of  density 
and  tracers  in  the  eddying  flows  cannot  be  given  by  only  the  Eulerian  mean  circulation, 
but  the  distribution  also  involves  the  eddy-induced  advection,  a  process  similar  to  the 
Stokes  drift.  Therefore,  the  starting  point  is  the  time-mean  density  equation  that  takes 
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into  account  the  eddy  fluxes  as  follows: 


V-Vo  +  V(V'o')  = —  (2) 

dz 

where  o  is  the  time -mean  density  anomaly,  primes  denote  the  perturbations  from  this 
mean  due  to  transient  eddies,  and  B  represents  the  vertical  density  flux  due  to  small-scale 
processes  and  air-sea  fluxes.  While  we  mostly  associate  the  eddy  fluxes  with  the 
mesoscale  variability,  it  is  convenient  to  include  in  it  all  transient  processes,  such  as  the 
effects  of  the  seasonal  cycle  and  variation  in  the  mixed-layer  depth. 

The  eddy  fluxes  in  Eq.  (2)  can  be  broken  down  into  two  distinct  components 
(Radko  and  Marshall  2006): 

_ „* 

-  The  adiabatic  advection  by  the  eddy-induced  velocity  of  the  residual  mean  theory  ( V  ), 

_ „* 

which  is  assumed  to  be  non-divergent  (V  -V  =  0  ); 


-  The  remaining  diabatic  component  that  can  be  written  without  loss  of  generality  as  the 
vertical  convergence  of  the  vertical  flux 


r  SB,d^ 


dz 


As  a  result  of  this  transformation,  Eq.  (2)  becomes: 

o  dB 

V  res  '  VCr  = - 

dz 


(3) 


where  the  residual  velocity  Vres  =  V  +  V  represents  the  advection  of  density  and  tracers 
by  both  mean  field  and  adiabatic  eddies,  and  B  =  B  +  Bedd  includes  the  diabatic  effects  of 
the  small-scale  mixing  and  eddies. 

Following  Marshall  and  Radko  (2003),  we  separately  discuss  the  dynamics  of  the 
thin,  vertically  homogeneous  mixed-layer  ( —hm  <  z  <  0 )  and  the  stratified  interior 

In  the  mixed-layer  the  residual  buoyancy  equation  reduces  to: 


d<jm  do 
— — +  v  — 1 

^  res 

OX 


dy 


dB_ 

dz 


(4) 


where  om  is  the  mixed-layer  density  anomaly.  Hence,  the  integration  of  Eq.  (4)  over  the 
depth  of  the  mixed-layer  results  in: 
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u^r+v^r  =  B>+B^.-B 

ox  oy 


z=-h„. 


(5) 


u  u 

with  U  =  J  uresdz ,  V  =  J  vresJz ,  B(l  being  the  air-sea  density  flux,  Bmle  a  contribution 


—h„ 


—hm 


from  the  diabatic  eddies  in  the  mixed-layer,  and  B  the  vertical  density  flux  due  to 

z=-hm 


diabatic  processes  immediately  below  the  mixed-layer. 

In  the  mixed-layer  buoyancy  equation  (5),  the  total  density  flux  on  the  right-hand 
side  is  clearly  dominated  by  the  direct  air-sea  forcing  B0 .  The  estimates  in  Radko  and 

Marshall  (2006),  in  the  context  of  the  ACC,  suggest  that  the  sea  surface  flux  B0  exceeds 
Bmle  by  an  order  of  magnitude;  the  diabatic  eddy  fluxes  play  an  even  lower  role  on 
planetary  scales  (Radko  2006)  that  are  the  primary  focus  of  this  study.  In  addition,  both 
small-scale  mixing  and  the  diabatic  eddy  effects  are  greatly  reduced  below  the  mixed- 
layer  and  therefore  Bn  »  B\  .  Hence,  Eq.(5)  may  be  approximated  by: 

I  z=-hm 


U^s-  +  V^s-  =  Bn 


dx 


dy 


(6) 


The  integration  of  the  residual  continuity  equation  over  the  depth  of  the  mixed- 
layer,  associated  with  the  boundary  condition  that  the  vertical  component  of  the  residual 
velocity  vanishes  at  the  sea  surface,  leads  to: 


dU  8V  _  - 

dx  dy  res  z=~/!,“ 


(7) 


For  convenience,  in  following  equations,  wres  h  will  be  written  as  wres . 


The  next  step  consists  of  introducing  a  convenient  measure  of  the  strength  of  the 
residual  flow  that  enters  (escapes)  the  mixed-layer  at  a  given  sea  surface  density  range. 
Consider  a  thin  strip  of  fluid  bounded  by  two  nearby  isopycnal  surfaces  (a  and  cr  + Act) 
as  indicated  by  the  shaded  area  in  the  schematic  in  Fig.  4.  The  integral  of  wresover  this 
area  ( S )  is  defined  by: 

G(<x)  =  JJw„,rf5-  (8) 
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Figure  4.  Schematic  diagram  illustrating  the  proposed  diagnostic  framework. 


The  area  of  integration  in  Eq.  (8)  extends  from  one  intersection  with  the  coastal 
boundary  to  another  if  the  buoyancy  contours  are  blocked  by  the  land  [see  the  schematic 
diagram  in  Fig.  4],  and  it  otherwise  may  include  contributions  from  several  ocean  basins. 
For  closed  surface  density  contours,  for  instance  in  the  re-entrant  Antarctic  Circumpolar 
Current,  the  integral  pertains  to  the  entire  area  between  contours  cr  and  a  +  Acr . 


The  goal  of  the  following  derivation  is  to  express  G  in  terms  of  observed 
quantities.  By  using  the  relation  between  wres  and  U  and  V  provided  in  Eq.  (7),  and 

Green’s  theorem,  the  area  integral  G(cr)  can  be  expressed  by: 


=  JJ  ^,dS  =  III  | G  Ajtt  =  j>(Udy  -  Vdx) 


(9) 


as 


where  dS  is  the  boundary  of  the  region  S  .  Note  that  the  direction  for  the  loop  integral  is 
counterclockwise. 

Assuming  no  flux  across  the  coastal  boundaries  for  zonally  blocked  flows  (or  the 
exact  cancellation  of  the  zonal  fluxes  in  case  of  a  re-entrant  flow),  the  integral  in  Eq.  (9) 
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may  be  simplified  to: 


G[cr^  =  ^(Udy -Vdx)-  J  ( Udy-Vdx ) 


(10) 


<T+A<7 


where  the  integration  along  the  isopycnals  is  in  the  positive  x-dircction. 

By  introducing  the  unit  vector  nonnal  to  the  buoyancy  contours  n  = 

(10)  can  be  conveniently  rewritten: 

G (cr)  =  J  U -n  dl  -  Jf7 -n  dl 

cr+Acr  cr 


I  Vo-,, 


Eq. 


(ID 


with  U  =  (U,V)  and  /  being  the  arc  length  along  the  mixed-layer  density  contours.1 


Eq.  (11)  is  readily  interpreted  as  follows:  the  area- integrated  vertical  residual  flux  at  the 
bottom  of  the  mixed-layer  in  a  steady  state  is  equal  to  the  sum  of  the  lateral  residual 
fluxes  into  this  area.  This  statement  is  consistent  with  a  simplified  mass  budget  in  the 
volume  defined  by  the  two  isopycnal  surfaces  within  the  mixed-layer,  as  depicted  by  the 
schematic  diagram  in  Fig.  5.  This  mass  budget  is  here  expressed  as  a  simple  volume 
budget  since  the  density  is  almost  constant  within  the  considered  volume  and  the 
contribution  of  evaporation/precipitation  at  the  sea  surface  to  the  volume  balance  is 
neglected. 

Taking  the  limit  A  cr  — »  0 ,  Eq.  (11)  becomes: 


G(a)  =  —  -Act 

(12) 

where 

F(cr )  =  Jt7  •  n  dl 

(13) 

cr 


1  While  the  schematic  in  Fig.  4  and  the  related  discussion  represents  the  case  in  which  density 
increases  northward  (the  configuration  relevant  for  the  Northern  hemisphere),  it  can  be  readily  shown  that 
our  formulation  applies  to  both  Northern  and  Southern  hemispheres. 
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Mixed-layer 

depth 


Figure  5.  Physical  interpretation  of  Eq.  (11):  mass  budget  in  the  volume  delimited  by  the 
isopycnal  surfaces  a  and  cr  +  Acr  within  the  mixed-layer. 


Eq.  (6)  may  be  written  as  U  n=  - — —  and  then  substituted  in  Eq.  (13)  to  obtain 

| 


the  following  expression  of  F : 

F(c r)  =  [  ii/ 
V<T 

<7  |  m 


(14) 


Since  the  air-sea  fluxes  and  the  sea  surface  density  are  easily  accessible  from  data,  the 
preceding  analysis  offers  a  simple  procedure  for  estimating  the  total  residual  volume  flux 
that  enters  (escapes)  the  mixed  layer  from  (into)  the  ocean  interior  within  the  surface  area 
S'  (crj  <  cr  <  cr2 )  that  is  bounded  by  the  two  density  contours  a1  and  a2  (cr2  >  oy) : 

<r2  Q(q-\ 

V(a1,a2)=  JJ  wresdS  =  J  — ^ —  da  =F{ct2)-F{ux)  (15) 

S[<jl<a<a2 )  (J\  ^ 


where  ^(cr)  is  given  by  Eq.  (14). 
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B.  ALTERNATIVE  METHOD 

Instead  of  computing  F(a) along  the  ispopycnals  cr,  and  cr2,  V(avcx2)may  be 

obtained  by  an  alternative  method  using  the  area  integration  over  the  surface  bounded  by 
two  isopycnals.  This  alternative  technique  is  similar  to  that  in  Walin’s  (1982) 
formulation. 


By  introducing  the  expression  for  F  [Eq.  (14)]  into  Eq.  (11),  we  evaluate  the 
vertical  residual  flow  that  enters  (escapes)  the  mixed-layer  over  the  area  S  (defined  above 
and  sketched  in  Fig.  4): 


G{a)  =  \\wresdS  =  F(a  +  \a)-F(a)  =  J 

S  ct+Ao-  | V  °m  a  | V 


B„ 


41 


(16) 


cr. 


Alternatively,  Eq.  (16)  may  be  expressed  as: 


G(o-)  =  —  f 

v  '  An  J 


Bn^-xH-jBnT^rdl 


\cr+A(j 


V  cr 


lV 


cr, 


m\  J 


=  ^-f  J  B0Ldl  - 1  B0Ldl 


<  cr+Acr 


(17) 


where  L  =  - - -  is,  approximately,  the  distance  between  the  two  isopycnal  contours  cr 

Vcr  ‘  ' 

I  m  I 

and  cr  +  Acr .  Therefore,  Eq.  (17)  becomes: 


G(a) 


1 

Act 


JJ  B0  (cr  +  A a)dS  -  JJ  B0  (a)dS 

s  s  J 


(18) 


Taking  the  limit  Acr  — »  0 ,  we  arrive  at: 


GM  =  JR .dS*4>  Jj  B„dS 


8az 


Act 


(19) 


and  finally  obtain  the  alternative  expression  fork (crj,cr2)  : 


^G(a)  8 

V(a1,a2)=  J  da  =  x 

J  A  G  8(7 

( 

rr 

8 

( 

rr 

JJ  B«ds 

da 

JJ  B°ds 

Cl 

J 

<7=<T2 

VA(o-) 

A 

(20) 


13 


where  5  (cr)  is  the  area  bounded  by  the  mixed-layer  density  contour  <jm  =  <j  and  a 
reference  isopycnal  a0  located  in  the  same  hemisphere  and  verifying  a0  <  a .  Note  that 
the  exact  value  of  that  reference  density  does  not  affect  Eq.  (21);  therefore,  for 
convenience,  we  may  define  a0  as  the  minimum  density  with  distinct  non-intersecting 
northern  and  southern  contours,  as  depicted  in  Fig.  6a. 

The  alternative  method  leads  to  the  following  expression  for  F: 


While  this  expression  is  identical  to  the  water  mass  transformation  function 
(Walin  1982),  its  role  and  interpretation  in  the  model  derived  here  are  less  conventional. 
First,  we  note  that  unlike  Walin’s  model  and  its  previous  extensions,  this  fonnulation 
takes  into  account  the  volume  transfer  driven  by  the  adiabatic  component  of  mesoscale 
eddies.  More  importantly,  in  this  study,  function  F{o )  is  defined  [Eq.  (13)]  in  terms  of  a 
horizontal  volume  flux  in  the  mixed-layer  crossing  a  given  isopycnal  surface  cr.  This 
connection  makes  it  possible  to  calculate  the  net  residual  volume  flux  that  escapes 
(enters)  the  mixed-layer  from  (into)  the  ocean  interior  within  an  area  bounded  by  two 
density  contours  [Eq.  (15)  and  Eq.  (20)].  Walin  (1982)  did  not  separate  the  water  mass 
transfonnation  within  the  mixed-layer  from  that  occurring  in  the  ocean  interior. 

It  is  also  essential  to  note  that  the  foregoing  model  makes  it  possible  to  distinguish 
the  density  forcing  in  the  Northern  and  Southern  Hemispheres.  On  the  contrary,  the 
original  Walin’s  formulation  involved  the  total  transformation  at  a  given  interior 
isopycnal  surface,  in  which  contributions  from  the  air-sea  density  fluxes  in  both 
hemispheres  were  grouped  together. 
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C.  DIAPYCNAL  VOLUME  FLUX 

As  defined  above,  a0  denotes  the  minimum  value  of  mixed-layer  density  for 

which  the  studied  basin  contains  two  non-intersecting  isopycnal  contours.  As  shown  in 
Fig.  6a,  for  every  cr  >  cr0  there  are  two  distinct  values  of  F  corresponding  to  the  southern 

and  northern  density  contours;  these  are  denoted  as  Fs{cr )  and  FN(a)  respectively. 


Figure  6.  Schematic  representation  of  the  density  structure  in  a  zonally  bounded  ocean. 

(a)  Mixed-layer  density  distribution.  For  each  value  of  cr  >  cr0  there  are  two  non¬ 
intersecting  contours  of  surface  density;  one  in  the  northern  part  (indicated  by  the 
gray  shading)  and  one  in  the  southern  part  (not  shaded),  (b)  The  vertical 
meridional  cross-section  of  density.  The  residual  flux  entering  the  density  layer 
bounded  by  the  isopycnals  cr,  and  cr2  from  the  mixed  layer  (SVs  +SVN)  is 

balanced  by  the  diapycnal  flux  [  Vd ( a2 )-Vd(crl)]. 


Since  the  residual  flow  is  assumed  to  be  non-divergent,  the  total  residual  flux 
(Vd)that  crosses  the  isopycnal  cr  in  the  ocean  interior  (see  the  schematic  in  Fig.  6b)  is 
exactly  equal  to  the  flux  from  the  ocean  interior  into  the  mixed-layer  within  the  area 
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defined  by  the  condition  <jm  <  cr  : 

vd(°)=  \\  G(c?')d(j' =  Fs(cr)  +  FN(<j)  fova>cr0.  (22) 

The  sign  convention  in  Fig.  6  is  such  that  Vd  is  positive  (negative)  for  an  upward 
(downward)  integrated  flux  across  the  isopycnal  in  the  ocean  interior. 


D.  ADIABATIC  ADVECTION 

The  next  step  (and  the  key  strength  of  the  proposed  technique)  is  to  quantity  the 
magnitude  and  pattern  of  the  isopycnal  pole-to-pole  MOC  component,  caused  by  the 
adiabatic  interior  advection,  from  the  surface  data.  For  this,  the  residual  volume  flux 
entering  (escaping)  the  mixed-layer  from  (into)  the  interior  is  computed  separately  in  the 
Northern  and  Southern  Hemispheres. 

The  highest  value  of  potential  density  that  can  be  found  in  the  mixed-layer  is 
denoted  by  crmax ,  and  cr  is  an  arbitrary  value  within  the  interval  [cr0,crmax],  For  each 
hemisphere,  Eq.  (15)  is  used  to  express,  in  tenns  of  F  ,  the  net  residual  flux  at  the  bottom 
of  the  mixed-layer  within  a  density  range  of  [cr,  crmax  1 ,  as  follows: 


< 


VN  (<?)  =  jj  WresdS  =  Fn  (^max  )  -  F„  (a) 

north 

vs  (<?)  =  JJ  wresds  =FS  (o-max  )  -  Fs  (a) 

S(a<a'<amax ) 
south 


(23) 


Since  the  residual  flow  is  non-divergent,  the  sum  of  all  fluxes  entering  the  interior 
volume  bounded  by  the  isopycnals  crmax  and  cr  is  zero  [see  the  schematic  in  Fig.  6b]: 


K  (omax  )-Vd(<j)-VN(<j)-Vs(<j)  =  0  (24) 

Note  that  in  the  absence  of  any  cross-isopycnal  flows,  Eq.  (21)  would  reduce  to 
Vn(<t)  =  -Vs(<t),  which  simply  means  that  all  the  water  subducted  in  the  Southern 
(Northern)  Hemisphere  upwells  in  the  Northern  (Southern)  one.  In  this  case,  the 
isopycnal  flux  throughout  the  density  layer  [cr,crmax]  is  unifonn  and  the  adiabatic  volume 

transport  is  F  (cr)  =  Vn(<j)  =  -Vs (cr)  . 
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However,  in  reality,  the  diapycnal  flux  across  any  isopycnal  surface  is  not  zero, 
and  consequently  VN(a )  is  not  equal  to  Vs(a )  .  Therefore,  the  isopycnal  flux  cannot  be 
precisely  detennined  from  Vs  and  VN ,  and  we  need  to  define  a  way  to  assess  the 
magnitude  of  Va.  It  is  reasonable  to  assume  that  representative  values  of  Va  lie  in 
between  VN(a )  and  - Vs(a )  .  Thus,  without  any  additional  assumptions,  the  simplest  and 
most  consistent  estimate  of  the  isopycnal  flux  Va  within  the  density  range  [cr,crmax]  is 

given  by  the  linear  average  of  the  two  net  residual  fluxes  at  the  bottom  of  the  mixed- 
layer: 

W=2[W-W]  <25> 

Note  that,  with  this  definition,  a  positive  value  of  Va  is  associated  with  a 

downward  residual  flux  at  the  base  of  the  mixed  layer  in  the  Southern  Hemisphere  and  its 
corresponding  upward  flux  in  the  Northern  Hemisphere.  While  this  definition  seems 
somewhat  arbitrary  at  first  sight,  its  application  to  the  global  overturning  circulation 
model,  described  in  section  IV-D,  indicates  that  it  is  surprisingly  accurate  in  predicting 
the  actual  inter-hemispheric  volume  exchange.2 

Since  the  length  of  the  mixed-layer  density  contours  inevitably  reduces  to  zero  as 
the  value  of  density  approaches  its  highest  possible  value  ( cr  — »  crmax ),  we  may  assume 

that  Fv(crmax)  «  fs(crmax)  «  0  in  Eq.  (23),  and  therefore  Eq.  (25)  reduces  to: 

KM  =  ^[Fs(«)-Fn(*)]  (26) 

It  is  important  to  note  that  the  expression  above  is  most  relevant  for  the  case  of  a 
global  ocean  circulation.  On  the  contrary,  for  computations  within  each  basin,  where  we 
only  consider  the  isopycnal  contours  zonally  bounded  by  land,  a  better  way  to  define  Vs 
and  VN  than  that  expressed  in  Eq.  (23)  may  be: 

2  Note  that  cr  in  the  preceding  formulation  cannot  be  less  than  cr0 ,  and  therefore  Eq.  (25)  does  not 
describe  the  overturning  in  the  upper  tropical  thermocline  where  cr  <  crfl .  However  the  region 
corresponding  to  this  density  range  occupies  only  a  small  fraction  of  the  thermocline,  with  a  maximum 
depth  of  200  m,  and  thus  its  relative  contribution  to  the  overall  overturning  is  expected  to  be  negligible. 
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VN{°)=  JJ  wn-dS  =  Fn  (°')~  Fn (an ) 

north 

Vs(a)=  JJ  WresdS  =FS  O0 

5’(<J0<cr'<cr) 
south 


With  this  definition,  the  integration  is  carried  out  over  the  area  defined  by  the 


mixed-layer  density  range  [cr0,cr],  which  is  more  representative  of  the  considered 


oceanic  basin  than  the  density  interval  [cr,crmax  ] .  Moreover,  with  the  alternative 


definition  of  Vs  and  VN  provided  in  Eq.  (27),  the  adiabatic  flux  Va  defined  in  Eq.  (25) 
now  represents  the  isopycnal  volume  transport  within  the  density  range  [cr0,cr] . 


E.  ADIABATIC  FLUX  RATIO 

Va  (cr)  in  Eq.  (26)  represents  the  integrated  isopycnal  volume  flux  over  the 
density  interval  [cr,crmax] .  Therefore,  the  fraction  of  this  adiabatic  volume  flux  within  the 

8V 

density  layer  \<j,<j  +  Acrl  is  — -  Act  ,  where  Act  is  the  density  increment. 

do 

One  of  the  objectives  of  this  study  is  to  compare  the  diapycnal  and  isopycnal 
components  of  the  overturning  circulation,  and  some  quantification  of  their  relative 
strengths  is  obtained  through  the  following  parameter: 


r  = 


(28) 


(SV.\ 

=  rms 

sva 

\do/ 

<J0<(T<(Jmax 

do 

rms 

<j0<cr<<Trn. 


Vd(<r)  I,  and  A  om  =oimx  -o„. 


Parameter  y  is  referred  to  hereafter  as  the  adiabatic  flux  ratio.  If  y  » 1 ,  the 
interior  flow  is  largely  adiabatic,  whereas  y  « 1  indicates  the  importance  of  the  interior 
diffusive  processes  for  maintenance  of  the  MOC. 
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F.  EKMAN  FLUX 

As  defined  above,  the  residual  flux  includes  the  mean  circulation  and  the  eddy- 
contribution.  The  F.kman  flux,  which  is  generated  by  the  mechanical  forcing  by  the 
surface  winds,  dominates  the  lateral  mean  flux  in  the  mixed-layer.  Therefore,  in  regions 
where  the  eddy  activity  is  weak,  the  residual  flux  may  be  almost  equal  to  the  Ekman  flux. 

In  particular,  the  Ekman  pumping  velocity  at  the  bottom  of  the  Ekman  layer  is 
given  by: 


d_ 

dx 


(  T  ^ 
E>  y 

d 

fT  3 

L0x 

=  curl 

(  -  \ 

z 

yPof  J 

8y 

\Pof  J 

\  Pof  V 

(29) 


where  r0v  and  r0  are  respectively  the  zonal  and  meridional  components  of  the  wind  stress 
f  at  the  sea  surface,  p0  is  the  reference  seawater  density,  and /is  the  Coriolis  parameter. 


In  a  similar  way  to  the  expression  of  the  net  residual  flux  at  the  bottom  of  the 
mixed-layer  [Eq.  (23)],  we  may  define  for  each  hemisphere  the  net  Ekman  flux  at  the 
bottom  of  the  Ekman  layer  within  a  density  range  of  [cr,  crmax  ] ,  as  follows: 


^N-EK  (°")  JJ  WEdS 

5(o-<o-'<CTmax) 

north 


< 

^ S-EK  (cr)  —  JJ  WEdS 

S(c<o-'<crmax) 

south 


(30) 


As  mentioned  in  Section  II-D,  in  the  case  of  a  specific  oceanic  basin,  where  only 
the  isopycnal  contours  zonally  bounded  by  the  land  are  considered,  we  estimate  the 
overturning  circulation  in  the  density  range  [cr0,cr].  Therefore,  the  alternative 
expressions  of  the  northern  and  southern  Ekman  fluxes  are  in  this  case: 

Vn-ek  (cr) —  JJ  wEdS 


< 


5,(<t0<<t'«t) 

north 

^.S'-EK  (Cr)  —  JJ  WEdS 


(31) 


5,(<t0«t'«t) 

south 
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Finally,  by  analogy  to  the  previous  definitions  relative  to  the  residual  flux,  we 
may  define  the  Ekman  adiabatic  advection  (wind-driven  component  of  Va): 

K-BK  M  \ [VN-EK  W  ■ - : VS-BB  W]  (32) 

The  comparison  of  Va_EK(cr)  and  Va(cr)  allows  us  to  assess  the  role  of  the  wind 
stress  in  the  establishment  and  maintenance  of  the  MOC. 
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III.  PARAMETER  DEFINITION  AND  DATABASE  DESCRIPTION 


The  technique  described  in  the  previous  chapter  offers  a  quantitative  description 
of  the  two  components  (adiabatic  and  diapycnal)  of  the  MOC  based  on  sea  surface 
observations.  The  next  step  consists  of  choosing  the  datasets  in  order  to  provide  an 
accurate  assessment  of  the  oceanic  overturning. 

Since  the  MOC  varies  on  time-scales  which  greatly  exceed  the  typical  scale  of  the 
air-sea  fluxes  (diurnal  cycles),  the  model  used  in  this  study  must  rely  on  long-term  means 
of  observed  parameters,  and  the  computations  carried  out  below  use  climatology  datasets. 
This  approach  is  consistent  with  the  variability  of  the  studied  phenomena  and  secondly  it 
allows  the  analysis  of  climate  variations,  which  is  one  of  the  main  motivations  for 
studying  the  MOC. 

The  computations  in  this  study  use  annual  means  for  all  quantities  and  it  is 
essential,  throughout  the  data  description,  to  examine  whether  the  climatology  averages 
accurately  represent  the  reality. 


A.  MIXED-LAYER  DENSITY 

1.  Definition 

One  of  the  quantities  used  in  the  computation  of  the  residual  flux  is  the  mixed- 
layer  density  anomaly  ( <jm  ).  The  mixed-layer  density  is  conceptually  well-defined:  the 

intense  mixing  within  the  mixed-layer  (ML)  makes  the  seawater  density  vertically 
uniform.  Therefore,  in  the  conceptual  models,  the  density  field  at  any  depth  above  the 
pycnocline  accurately  represents  the  ML  density  distribution. 

However,  in  reality,  the  annual-mean  density  is  not  vertically  uniform  in  the 
mixed-layer  and  the  choice  of  one  particular  depth  to  characterize  the  ML  properties  is 
not  trivial.  This  choice  is  made  more  difficult  by  the  large  seasonal  and  spatial  variability 
of  the  ML  depth.  Indeed,  our  theory  uses  a  spatially  uniform  model  of  the  mixed-layer 
depth  for  the  computation  of  the  density  gradient  that  is  needed  in  the  calculation  of  the 
MOC  components  [Eq.  (14)]. 
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Fig.  7  shows  the  spatial  distribution  of  the  ML  depth:  the  contours  represent  the 
annual  means  provided  by  the  Levitus  94  climatology.  This  map  indicates  the  large 
difference  between  the  subtropical  and  polar  regions.  These  differences  are  significantly 
reduced  here  because  the  plotted  data  are  the  annual  averages  and  this  filters  out  the 
seasonal  variations.  The  observations  in  Fig.  7  suggest  that  the  mixed-layer  density  can 
be  approximated  by  the  density  at  the  sea  surface.  This  idea  is  supported  by  the  plot  in 
Fig.  8,  which  gives  the  annually-  and  zonally-averaged  ML  depth  with  its  temporal 
standard  deviation  computed  over  the  monthly  mean  values. 

On  the  other  hand,  as  shown  in  Table  1,  the  temporal  variations  of  density  tend  to 
significantly  increase  as  the  chosen  depth  is  closer  to  the  sea  surface.  Thus,  in  order  for 
a(h)  to  accurately  represent  the  mixed-layer  density  distribution  throughout  the  year, 

we  have  to  choose  depth  h  large  enough  to  have  a  low  temporal  standard  deviation. 
Moreover,  a  deeper  level  within  the  ML  is  more  relevant  for  describing  the  density  layer 
of  the  MOC  in  the  ocean  interior. 

Thus,  a  compromise  must  be  found  in  order  to  select  a  depth  which  is  almost 
always  and  everywhere  within  the  mixed-layer  and  at  which  the  density  temporal 
variability  is  relatively  low.  The  depth  chosen  for  this  study  is  40  meters,  based  on  Fig.  8 
and  Table  1.  Note  also  that  the  40  m  depth  is  not  within  the  ML  only  in  the  equatorial 
region  and  in  the  Arctic  Ocean.  The  equatorial  area  is  excluded  from  our  analysis:  we 
only  consider  regions  where  cr  >  cr0  and  the  technique  derived  in  the  previous  chapter 

does  not  describe  the  overturning  in  the  upper  tropical  thennocline  where  a<  a0.  Since 

the  North  Pole  region  is  mostly  covered  by  ice,  the  density  forcing  at  the  sea  surface  there 
will  not  be  taken  into  account  in  our  computations. 
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Figure  7.  Annual  mean  of  the  oceanic  mixed-layer  depth  given  in  meters  and  based  on 
Levitus  94  climatology.  White  areas  represent  the  continents  and  coastal  regions 
where  the  mixed-layer  depth  is  less  than  20  m. 


In  order  to  more  accurately  represent  the  ML  density,  the  small  effect  of  the 
pressure  on  the  density  must  be  removed;  hence,  the  potential  density  anomaly  at  40  m  is 
used  for  the  computations.  The  potential  density  anomaly  (cr0)  is  related  to  the  potential 

density  ( pg  in  kg  m  3 )  as  follows: 

Pg  =  1000  +  <Jg 

Note  that  using  the  potential  density  field  to  represent  the  ML  density  brings  our  analysis 
closer  to  the  conceptual  models  of  the  mixed-layer,  since  this  parameter  has  a  lower 
vertical  gradient  than  density. 
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Figure  8.  Annual  mean  and  temporal  standard  deviation  of  the  zonally-averaged  mixed- 
layer  depth.  The  plotted  data  are  based  on  Levitus  94  climatology.  The  solid  line 
indicates  the  zonally-  and  annually-averaged  ML  depth  distribution  and  the  dotted 
lines  point  out  its  temporal  standard  deviation  computed  from  the  monthly  ML 
depth  means.  The  dashed  horizontal  line  highlights  the  40  m  depth. 


Depth  (m) 

Standard  deviation  of  potential  density 
(kg  m'3) 

0  (surface) 

4.7945 

10 

0.3995 

20 

0.3565 

30 

0.3184 

40 

0.2670 

50 

0.2403 

75 

0.1720 

100 

0.1349 

Table  1.  Mean  temporal  standard  deviation  of  potential  density  at  each  characteristic  depth 
of  the  mixed-layer.  These  statistical  values  are  computed  from  the  Levitus  94 
monthly  means  of  potential  density  at  all  grid  points  of  the  whole  ocean. 
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2.  Levitus  Climatology 

The  oceanic  database  used  in  this  study,  for  the  potential  density  distribution  at  40 
m,  as  well  as  temperature  and  salinity  needed  in  the  calculations  of  the  density  flux  at  the 
sea  surface,  is  Levitus  94  climatology.  This  database  consists  of  a  world  atlas  of 
objectively  analyzed  mean  fields  of  major  oceanic  parameters  at  the  annual,  seasonal,  and 
monthly  time  scales. 

All  the  analyzed  fields  are  offered  on  a  one-degree  latitude-longitude  grid  at 
selected  standard  depth  levels  from  the  sea  surface  to  5500  m  depth.  The  standard  levels 
relative  to  the  ML  are  0,  10,  20,  30,  50,  75,  100,  125,  150,  200,  250,  and  300  meters. 
Therefore,  for  the  40  m  potential  density  anomaly  distribution,  we  have  linearly 
interpolated  the  30  m  and  50  m  analyzed  fields. 

The  atlas  datasets  have  been  produced  by  the  Ocean  Climate  Laboratory  (OCL)  at 
the  National  Oceanographic  Data  Center  (NODC),  which  is  supported  by  the  National 
Oceanic  and  Atmospheric  Administration  (NOAA)  Climate  and  Global  Change  program 
to  provide  quality  controlled  oceanographic  databases.  In  particular,  they  include  an 
objective  analysis  of  all  the  in  situ  measurements  for  the  period  1900-1992. 

The  Levitus  climatology  is  considered  by  all  oceanographers  to  be  one  of  the 
most  reliable  databases  for  the  oceanic  parameters.  It  is  accessible  through  the  website 
http://iridl.ldeo.columbia.edU/SOURCES/.LEVITUS94  in  NetCDF  (Network  Common 
Data  Form)  format. 

3.  Potential  Density  Field 

The  annual  mean  distribution  of  potential  density  anomaly  at  40  meters 
(< <jg  (40m))  from  Levitus  94  climatology  is  given  in  Fig.  9.  The  isopycnal  contours  cr0  in 
both  hemispheres  are  highlighted  on  the  map.  Based  on  the  definition  of  cr0  given  in  the 
previous  chapter  and  the  monthly  mean  fields  of  <jg  (40 m)  [not  shown  here],  the  value  of 
cr0  has  been  estimated  at  a0  =  24.2 .  This  value  corresponds  to  the  minimum  potential 
density  with  distinct  non-intersecting  northern  and  southern  contours  over  the  whole 
ocean  and  throughout  the  year  (i.e.,  on  every  chart  of  <jg  (40 in)  monthly  means). 
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Figure  9.  Annual  mean  distribution  of  potential  density  anomaly  at  40  m  depth.  The 
dimensionless  displayed  data  have  been  computed  from  Levitus  94  climatology. 
The  thick  solid  lines  point  out  the  northern  and  southern  contours  of  cr0 

(cr0=  24.2).  Dotted  lines  indicate  partitioning  of  the  world  ocean  into  the 
individual  basins. 


B. 


DENSITY  FLUX  AT  THE  SEA  SURFACE 


1.  Definition 

The  second  surface  quantity  used  in  the  computation  of  the  residual  flux  is  the  air- 
sea  density  flux  ( B  in  kg  m1s~x ).  As  shown  previously  in  Eq.  (1),  this  flux  is  composed 
of  two  terms,  the  thermal  part  (f?, )  and  the  freshwater  component)  f>2 ) : 


B ,  =- 


ccH 


1  -S 


C, 


b2  =  M 


(. E-P)S 
1  -s 


(33) 

(34) 
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Note  that  many  studies  dealing  with  fluxes  at  the  sea  surface  refer  to  the  buoyancy 
flux:  this  latter  may  be  obtained  from  the  density  flux  by  simply  multiplying  B  by  the 

constant  ,  with  g  being  the  gravity  constant  (9.8  ml  s2 ) and  p(]  the  standard 

Po 

seawater  density  (l025  kg  in  !  ) .  Indeed,  buoyancy  b  is  defined  as  b  =  — —  p ,  where  p 

Po 

is  the  seawater  density. 

The  parameters  Cp ,  a  and  [1 ,  which  are,  respectively,  the  specific  heat  capacity, 

the  thermal  and  saline  expansion/contraction  coefficients  of  seawater,  are  not  constant 
and  depend  on  temperature,  salinity,  and  pressure: 

a  =  a(T,S,p ) 

• P  =  P(T,S,p ) 

Cp=Cr(T,S,p) 

Therefore,  the  temperature  ( T)  and  salinity  (S)  distributions  at  the  sea  surface 
(where  the  hydrostatic  pressure  p  is  equal  to  0),  which  are  available  in  Levitus  94 
climatology,  will  be  used  for  determining  the  air-ocean  density  flux.  The  three 
coefficients  are  then  calculated  on  the  basis  of  the  UNESCO  algorithms  for  computation 
of  fundamental  properties  of  seawater  (1983). 

The  heat  flux  from  the  atmosphere  into  the  ocean  (If),  the  evaporation  (E)  and 
precipitation  ( P )  rates  are  obtainable  in  meteorological  databases.  In  particular,  the 
following  fields  may  be  found  in  the  different  climatological  databases: 

Net  solar  energy  flux  at  the  sea  surface  [Qs ) ; 

Net  longwave  radiation  at  the  sea  surface  [QL ) ; 

Latent  heat  flux  ( HL ) ; 

Sensible  heat  flux  [Hs ) ; 

Total  precipitation  at  the  surface (f5),  which  is  sometimes  divided  into  two 
fields  (large-scale  and  convective  precipitation). 

Note  that,  for  the  four  thermal  fluxes  (QS,QL,HL,HS) ,  the  convention  used  here 
associates  positive  values  to  upward  fluxes  (from  the  ocean  into  the  atmosphere). 
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The  five  parameters  described  above  allow  us  to  compute  the  heat  flux  from  the 
atmosphere  into  the  ocean  (H  in  W  m  1 )  and  the  evaporation  and  precipitation  difference 
{E-P  in  m  s  1 )  field  through  the  following  equations: 

H  =  -Qs-Ql-Hl-Hs  (35) 

and  E-P  =  — ^ — P  (36) 

PoLe 

where  Le  is  the  specific  latent  heat  for  evaporation  ( Le  «  2.5xl06  J  kg  1  j . 

2.  Atmospheric  Databases 

Four  different  databases  have  been  used  to  compute  H  and  E-P.  We  compare  their 
climatology  and  the  corresponding  predictions  of  the  MOC  components.  This  comparison 
answers  two  essential  questions:  the  quality  control  of  the  results  and  their  sensitivity  to 
the  observational  fields. 

•  The  European  Center  for  Medium-range  Weather  Forecast  (ECMWF)  Re- 
Analysis  (ERA-40)  is  a  global  analysis  of  many  observations  and  satellite  data 
streams  for  the  period  September  1957  through  August  2002.  This  dataset 
contains  56  surface  and  single-level  variables  on  an  equally-spaced  global  2.5 
degree  latitude-longitude  grid,  with  all  variables  reported  four  times  a  day.  In 
particular,  the  ERA-40  atlas,  the  web-version  of  which  was  released  on  6 
February  2006,  describes  the  climate  from  1979  through  2001,  which  is  the  period 
with  the  best  and  most  time-consistent  product  quality  for  the  entire  globe.  For 
this  reason,  this  study  uses  monthly  means  of  each  parameter  listed  above  from 
January  1979  through  December  2001.  The  data  may  be  downloaded  from 
http://dss.ucar.edu/datasets/dsl  18.0.  Note  that  the  data  are  provided  there  in  GRIB 
format  (binary  code  for  exchange  of  processed  data). 

•  The  joint  National  Centers  for  Environmental  Prediction  (NCEP)  and 
National  Center  for  Atmospheric  Research  (NCAR)  re-analysis  project  has 
performed  data  assimilation  using  historical  data  from  January  1948  to  the 
present.  Particularly,  the  NCEP-NCAR  Climate  Data  Assimilation  System  1 
(NCEP -NCAR  CDAS-1)  provides  monthly  means  of  each  needed  parameter  from 
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January  1949  up  to  the  current  month  on  an  equally-spaced  global  latitude- 
longitude  grid  containing  94x192  points.  In  order  to  be  consistent  with  the 
ECMWF  climatology,  we  use  the  subsets  from  January  1979  through  December 
2001  in  this  study.  These  datasets  may  be  obtained  at 
http://iridl.ldeo.columbia.edU/SOURCES/.NOAA/.NCEP-NCAR/.CDAS-l  in 

NetCDF  format  (Network  Common  Data  Form). 

•  The  Global  Air-Sea  heat  and  momentum  flux  Climatology  from  the 
Southampton  Oceanography  Centre  (SOC-GASC97)  provides  climatologic 
sea  surface  fluxes  based  on  in  situ  reports  within  the  COADS  1A  (Comprehensive 
Ocean  Atmosphere  Dataset  1A),  a  global  dataset  containing  of  order  30  million 
surface  observations  from  ships  an  buoys  collected  over  the  period  1980-1993.  It 
is  important  to  note  that  the  quality  of  the  fields  is  known  for  having  a  strong 
spatial  dependence  that  reflects  the  global  distribution  of  ship  observations. 
Quality  is  likely  to  be  high  in  the  well-sampled  North  Atlantic  and  North  Pacific 
oceans  but,  it  decreases  in  the  Southern  Hemisphere.  In  particular,  south  of  40  S, 
the  errors  in  the  fields  are  likely  to  be  large  and  the  existence  of  spurious  features, 
which  have  been  generated  during  the  objective  analysis  of  the  original  raw 
distributions,  must  be  pointed  out.  The  data  are  available  on  an  equally-spaced 
global  one-degree  latitude-longitude  grid,  in  NetCDF  format,  at 
http://iridl.ldeo.columbia.edU/SOURCES/.SOC/.GASC97.  The  dataset 
documentation  is  accessible  from  this  website. 

•  The  Da  Silva  atlas  of  Surface  Marine  Data  1994  (SMD94)  provides  the 
monthly  climatology  of  the  near-surface  atmospheric  properties  over  the  global 
ocean.  Covering  the  period  of  January  1945  through  December  1989,  this 
database  is  made  from  individual  observations  in  COADS  1A  and  an  objective 
analysis  using  Barnes’  weight  function,  which  is  the  same  response  function  used 
in  Fevitus  94  Climatology.  It  was  produced  at  the  University  of  Wisconsin- 
Milwaukee  (UWM)  by  A.  M.  da  Silva  and  C.  C.  Young  in  collaboration  with  S. 
Fevitus  from  the  National  Oceanic  and  Atmospheric  Administration  (NOAA)  / 
National  Oceanographic  Data  Center  (NODC).  The  datasets  are  available  on  an 
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equally-spaced  global  one-degree  latitude-longitude  grid,  in  NetCDF  format,  at 
http://iridl.ldeo.columbia.edU/SOURCES/.DASILVA/.SMD94.  The  dataset 
documentation  is  accessible  from  this  website. 


3.  Distribution  of  the  Density  Flux  at  the  Sea  Surface 

Based  on  Eqs.  (1),  (33),  (34),  (35),  and  (36),  the  density  flux  at  the  sea  surface  has 
been  computed  from  each  of  the  four  databases.  Figs.  10,  11,  12,  and  13  show  the  density 
flux  distributions  respectively  from  ECMWF  re-analysis,  NCEP/NCAR  re-analysis, 
SOC-GASC97,  and  DaSilva-SMD94  climatology.  In  each  figure,  the  thermal  and 
freshwater  components  (respectively  Bi  and  B2 )  are  also  displayed. 


Over  most  of  the  oceans,  the  density  flux  seems  to  be  driven  by  thermal  forcing. 
Figs.  10,  11,  12,  and  13  show  that  the  distributions  of  total  density  flux  B  and  its  thermal 
part  B]  look  similar.  Moreover,  by  looking  at  their  magnitudes,  we  see  that  B]  is 

generally  the  dominant  part  of B .  This  observation  is  supported  by  the  zonal  means  of 
annually-averaged  values  that  are  given,  respectively  for  each  meteorological  database,  in 
Figs.  14,  15,  16,  and  17. 


Thus,  from  a  global  prospective,  thermal  forcing  dominates  the  freshwater  forcing 
and  controls  the  density  flux  at  the  sea  surface.  However,  as  shown  in  Figs.  18,  19,  20, 


and  2 1 ,  which  provide  the  annually-averaged  distribution  of  log 
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for  each  database. 


the  zonal  variability  is  large.  In  some  places,  especially  in  subtropical  and  polar  regions, 
B2  greatly  exceeds  Bl .  These  observations  become  more  critical  when  we  consider 
seasonal  fields  [not  shown  here]. 


Figs.  10,  1 1,  12,  and  13  show  the  same  regions  of  high  density  flux:  these  are  in 
general  the  Western  Boundary  Currents  (WBC)  and  the  equatorial  band.  Indeed,  in  the 
regions  of  the  warm  Gulf  Stream,  Kuroshio  Current,  Agulhas  Current,  and  East 
Australian  Current,  the  heat  loss  by  the  ocean  into  the  atmosphere  is  large,  which  causes 
the  upper  ocean  to  gain  density.  These  currents  advect  heat  from  the  tropics  towards  the 
poles,  and  energy  exchanges  with  the  atmosphere  occur  in  their  paths. 
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Figure  10.  Sea  surface  density  flux  from  ECMWF  climatology.  Positive  values  correspond  to 
density  gain  by  the  ocean.  The  displayed  values  are  in  kg  in  2s  1  (a)  Total  density 
flux  B  (b)  Thermal  component  Bx  (c)  Freshwater  component  . 
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b) 


c) 


Figure  11.  Same  as  Fig.  10,  but  from  NCEP/NCAR  re-analysis  (same  period  as  ECMWF 
climatology:  1979-2001). 
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Figure  12.  Same  as  Fig.  10,  but  from  SOC-GASC97  climatology  (period:  1980-1993). 
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a) 


b) 


c) 


Figure  13.  Same  as  Fig.  10,  but  from  DaSilva-SMD94  climatology  (period:  1945-1989). 
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Along  the  equator,  the  upwelling  induced  by  the  Ekman  divergence  brings  cold 
water  to  the  upper  ocean  and  leads  to  vigorous  downward  heat  fluxes  from  the  warm 
tropical  air  above  the  surface.  The  negative  oceanic  density  flux,  caused  by  thennal 
forcing,  is  enhanced  by  the  large  precipitation  within  the  Inter-Tropical  Convergence 
Zone  (ITCZ)  band,  which  is,  on  average,  centered  on  the  equator. 

Concerning  the  freshwater  contribution,  the  large  evaporation  and  the  weak 
precipitation  in  subtropical  regions  generates  a  strongly  positive  contribution  to  the 
density  flux  at  the  surface.  These  regions  are  characterized  by  high  sea-level  pressure  and 
persistent  coverage  of  low-level  non-precipitating  clouds.  On  the  other  hand,  over  the 
equatorial  band  (mean  position  of  the  ITCZ)  and  in  the  sub-polar  areas,  precipitation  is 
greater  than  evaporation,  resulting  in  negative  oceanic  density  flux. 

The  four  databases  considered  in  this  study  provide  the  same  overall  picture  of  the 
surface  density  flux  distribution.  However,  significant  quantitative  differences  are  noted 
in  some  regions  for  the  density  flux  components  Bx  and  B2 . 

It  is  generally  agreed  that  ECMWF  and  NCEP/NCAR  are  the  two  most  reliable 
among  these  four  climatologic  databases.  A  number  of  studies  have  compared  these  two 
re-analysis  datasets,  and  many  of  them  (for  instance,  Bromwich  and  Fogt  2004)  have 
suggested  that  ERA-40  seems  to  have  better  skills  than  NCEP/NCAR  regarding  analyzed 
fields  after  1978  (the  beginning  of  the  modern  satellite  era).  This  may  be  explained  by  a 
better  assimilation  of  satellite  data  within  ECMWF  than  in  NCEP/NCAR,  which  is 
largely  constrained  by  the  station  observational  network.  Thus,  ECMWF  is  particularly 
useful  for  the  high  latitudes. 
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Figure  14.  Zonal  means  of  annually-averaged  density  flux  B  at  the  sea  surface  and  of  its  two 
components  Bx  and  B2  fromECMWF  climatology  (1979-2001). 


Latitude  (°) 

Figure  15.  Same  as  Fig.  14,  but  from  NCEP/NCAR  re-analysis  (1979-2001). 
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Figure  17. 
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Same  as  Fig.  14,  but  from  SOC-GASC97  climatology  (period:  1980-1993). 
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Same  as  Fig.  14,  but  from  DaSilva-SMD94  climatology  (period:  1945-1989). 
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Based  on  the  foregoing  considerations,  ECMWF  re-analysis  has  been  chosen  as 
the  reference  meteorological  database  in  this  study.  This  choice  is  supported  by  the  charts 
in  Fig.  22,  which  display  the  differences  between  the  annually-averaged  density  fluxes 
obtained  from  ECMWF  re-analysis  and  those  computed  from  the  other  datasets.  Indeed, 
ECMWF  appears  to  be  a  good  compromise  between  NCEP/NCAR,  DaSilva-SMD94,  and 
SOC-GASC97.  NCEP/NCAR  provides  higher  density  flux  than  ECMWF  by  the 
subtropical  west  coast  of  major  continents  and  over  a  large  part  of  the  subtropical  North 
Atlantic  Ocean.  On  the  other  hand,  DaSilva-SMD94  and  SOC-GASC97  give  smaller 
values  than  ECMWF  over  all  tropical  and  subtropical  oceans  expect  for  some  west  coast 
regions.  These  features  are  observed  for  almost  all  latitudes,  as  indicated  by  the  zonal 
means  of  these  density  flux  differences  in  Fig.  23. 

Table  2  provides  the  global  mean  value  of  the  density  flux  differences  to  assess 
the  global  accuracy  of  each  meteorological  database  versus  ECMWF.  These  values 
support  the  choice  of  ECMWF  as  the  reference  climatologic  database  for  this  study. 


Mean  density  flux  difference 
(kg  mV) 

B  ( NCEP  /  NCAR  )~B(  ECMWF ) 

1.8716xlCT7 

B(SOC  -  GASC97)  -  B  (ECMWF) 

-1. 3075x1 0”6 

B  ( DaSilva )  -  B  (ECMWF) 

-1.3756xl0~6 

Table  2.  Global  averages  of  density  flux  differences  between  ECMWF  re-analysis  and  the 
other  databases. 
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Figure  19.  Same  as  Fig.  18,  but  from  NCEP/NCAR  re-analysis  (1979-2001). 
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Figure  20.  Same  as  Fig.  18,  but  from  SOC-GASC97  climatology  (period:  1980-1993). 


Figure  21.  Same  as  Fig.  18,  but  from  DaSilva-SMD94  climatology  (period:  1945-1989). 
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Figure  22.  Differences  in  the  surface  density  flux  (a)  ^B[NCEP I NCAR^- B^ECMWF)^ 

(b)  [B(SOC-GASC97)-B(ECMWF)]  (c)  [b (DaSilva) -B(ECMWF)]. 
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Figure  23.  Zonal  means  of  the  differences  in  the  surface  density  flux: 

[B(NCEP I NCAR)- B{ECMWF)\,  [ B ( SOC - GASC97) -  B ( ECMWF)] ,  and 

[B (DaSilva) -B(ECMWF)] . 


C.  MOMENTUM  FLUX  AT  THE  SEA  SURFACE 

As  described  in  section  II-F,  the  computation  of  the  Ekman  flux  requires  the  zonal 
and  meridional  components  of  the  wind  stress  f  at  the  sea  surface,  denoted  by  r0t.  and 

r0  respectively.  Three  databases  have  been  used  for  the  wind  stress  components: 

ECMWF  re-analysis,  NCEP/NCAR  re-analysis,  and  DaSilva-SMD94  climatology,  with 
the  same  periods  of  record  as  those  defined  for  the  density  flux  calculations.  The 
computed  Ekman  flux  provides  a  convenient  check  of  the  residual  flux  in  regions  where 
the  eddy  activity  is  weak. 
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IV.  RESULTS  AND  ANALYSIS 


A.  GLOBAL  OCEAN 

Fig.  24  plots  the  zonally-averaged  meridional  cross-section,  in  which  the 
climatological  potential  density  anomaly  contours,  computed  from  Levitus94  database, 
are  displayed.  This  cross-section  is  similar  to  the  meridional  cross-sections  of  individual 
oceans  and  shows  a  density  distribution  that  is  analogous  to  the  schematic  in  Fig.  6b. 


Figure  24.  Zonally-averaged  potential  density  anomaly  (dimensionless)  for  the  upper  1000 
meters  of  the  ocean.  The  data  have  been  computed  from  the  Levitus  94 
climatology. 


An  advantage  of  applying  the  diagnostic  model  to  the  global  configuration  is  that, 
unlike  the  individual  basins,  the  water  mass  budget  of  the  global  ocean  does  not  contain 
the  contribution  from  the  external  sources  or  sinks.  The  density  range  considered  for  the 
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full  ocean-wide  residual  flux  budget  extends  from  cr0  =  24.2 ,  determined  in  the  previous 


chapter,  to  crmax  =  28.3 ,  the  highest  density  of  the  ice-free  sea  surface  area. 


1.  Diapycnal  Volume  Flux 

The  diapycnal  interior  flux  ( Vd ),  computed  from  the  four  meteorological 
databases  using  Eq.  (21)  and  Eq.  (22),  is  plotted  in  Fig.  25.  Its  general  pattern  and 
average  value  of  3.5  Sv  (one  Sverdrup  is  defined  by  1  Sv  =  106  in  s  1 ),  obtained  from  the 
ECMWF  climatology  (reference  database),  are  consistent  with  the  earlier  estimates  by  the 
Walin-type  calculations  (Donners  et  al.  2005). 

Forcr<26.1,  Vd(a )  is  positive,  which  implies  an  upward  integrated  flux  across 

the  cr  isopycnal  surfaces  in  the  ocean  interior.  This  result  is  a  consequence  of  the  strong 
low-latitude  upwelling  and  the  net  upward  air-sea  heat  flux  in  mid-latitude  regions. 
Indeed,  as  indicated  in  Eq.  (14),  a  positive  density  flux  into  the  ocean  (induced  by  a  net 
upward  air-sea  flux)  leads  to  a  positive  value  ofF  (cr)  in  each  hemisphere,  which  implies 
water  mass  transformation  to  higher  densities  within  the  mixed-layer. 

The  diapycnal  flux  Vd  then  changes  sign  and  remains  negative  over  the  density 
anomaly  interval  [26.1,  27.6].  This  may  be  interpreted  as  a  convergence  of  the  residual 
flow  within  the  mixed-layer  in  this  density  range.  The  negative  values  of  Vd  are 

associated  with  the  strong  sinking  motion  of  Antarctic  Intermediate  Waters  (AAIW)  and 
Sub-Antarctic  Mode  Waters  (SAMW),  which  occurs  along  the  Antarctic  Circumpolar 
Current  (ACC),  characterized  by  a  potential  density  anomaly  of  about  27;  indeed,  the 
cross-isopycnal  flux  Vd  has  a  negative  peak  value  at  cr  «  27. 1 . 

For  cr  >  27. 1 ,  Vd  increases  with  cr  and  becomes  positive  from  cr  =  27.6  .  In  fact, 
the  strong  downwelling  along  the  ACC  is  compensated  by  the  large-scale  upwelling  of 
North  Atlantic  Deep  Water  (NADW)  and  North  Pacific  Deep  Water  (NPDW),  which 
takes  places  in  the  Southern  Ocean,  between  the  ACC  and  the  Antarctic  continent.  Thus, 
our  diagnostic  results  in  Fig.25  are  consistent  with  the  distribution  of  water  masses  over 
the  whole  ocean,  sketched  in  Fig.  2. 
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Figure  25.  Diapycnal  interior  flux  Vd  in  the  global  ocean  computed  from  the  four 
meteorological  databases. 


Fig.  25  also  highlights  the  strong  discrepancies  among  the  Rvalues  obtained 
from  the  four  climatologic  databases.  Only  ECMWF  and  NCEP/NCAR  datasets  result  in 
Vd  profiles  consistent  with  the  previous  results  (Donners  et  al.  2005).  Flowever, 

NCEP/NCAR  values  are  comparatively  too  high  at  low  density  anomalies  (cr  <  26) . 

SOC-GASC97  and  Da  Silva  SMD94  databases  generally  provide  inconsistent 
profiles.  This  result  could  be  expected,  since,  as  described  in  the  previous  chapter,  the 
air-sea  fluxes  in  these  two  databases  are  questionable  for  a  large  part  of  the  whole  ocean. 

Table  3  gives  the  average  cross-isopycnal  flow  differences  between  values 
obtained  from  ECMWF  re-analysis  and  those  computed  from  the  other  databases,  over 
the  potential  density  anomaly  interval  [24.2,  28.3].  This  table  should  be  compared  with 
Table  2,  which  provides  the  global  averages  of  density  flux  differences  between  ECMWF 
climatology  and  the  other  databases.  The  comparison  shows  that  the  diapycnal  volume 
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flux  estimates  are  critically  dependent  on  the  density  flux  accuracy.  Indeed,  as  shown  in 
Fig.  26,  the  absolute  diapycnal  volume  flux  accuracy  is  empirically  related  to  the  absolute 
mean  density  flux  accuracy,  by  the  following  equation: 

VDdiff  =  5.38xl025  Bdifft  -7.48xl019  Bdiff 2  +3. 90x1 013  \Bdiff\  (37) 

where  VDdiff  is  the  average  diapycnal  volume  flux  accuracy,  Bdiff  is  the  mean  density 

flux  accuracy,  and  all  the  coefficients  in  Eq.  (37)  are  measured  in  the  corresponding  MKS 
units. 


Mean  diapycnal  volume  flux 
difference  (mV1) 

VD  (NCEP INCAR)  -  VD  (ECMWF) 

5. 0470x1 06  (  »  5  Sv) 

VD  ( SOC  -  GASC97)  -  VD  (ECMWF) 

-4.3530xl07  (*-43.5  Sv) 

VD  (DaSilva)  -  VD  (ECMWF) 

-5.2320xl07  (*-52  Sv) 

Table  3.  Means  of  diapycnal  volume  flux  differences  between  ECMWF  re-analysis  and  the 
other  databases. 
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Figure  26.  Empirical  relation  between  the  absolute  mean  density  flux  accuracy  and  the 
absolute  average  diapycnal  volume  flux  accuracy.  The  red  dots  represent  the 
differences  in  the  density  flux  means  and  in  the  associated  diapycnal  volume  flux 
averages  between,  from  left  to  right,  NCEP/NCAR  and  ECMWF,  SOC-GASC97 
and  ECMWF,  and  Da  Silva  SMD94  and  ECMWF,  respectively. 


2.  Adiabatic  Advection 

Fig.  27  gives  the  distributions  of  the  along-isopycnal  interior  flux  ( Va ),  computed 
from  the  four  meteorological  databases  using  Eqs.  (21),  (23),  and  (25).  The  plots  display 
the  residual  (total)  fluxes  and  their  Ekman  components.  These  plots  should  be  interpreted 
with  caution,  since  they  represent  the  integrated  isopycnal  interior  flux  over  the  potential 
density  anomaly  range  [cr,  crmax  ] .  Thus,  changes  in  the  direction  of  the  MOC  circulation 

(northward  or  southward)  are  observable  by  maxima  and  minima  in  the  Va  profiles. 


The  direction  of  the  overturning  at  each  cr  is  detennined  by  sign  of 


da 


Since  Va  is  defined  by  Va(a )  =  ^\_VN(a)-Vs(a)]  [Eq.  (25)],  an  increasing  (decreasing) 


Va  from  crmax  down  to  a0  corresponds  to  a  northward  (southward)  flux. 
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Therefore,  the  total  Va  pattern  makes  it  simple  to  identify  the  individual 
components  of  the  isopycnal  MOC.  Most  notably,  the  southward  flux  of  the  North 
Atlantic  Deep  Water  (NADW)  is  reflected  in  the  strong  decrease  in  Va  from  about 

<7  =  27.6  to  cr  =  27.1.  The  northward  moving  water  masses  originating  in  the  ACC 
(Antarctic  Intennediate  Water  and  the  Sub-Antarctic  Mode  Water)  also  have  a  clear 
signature  through  the  significant  increase  in  Va  from  about  a  =  27. 1  to  a  =  26. 1 . 

Based  on  the  residual  (total)  isopycnal  flux  profile  computed  from  ECMWF 
climatology,  the  sketch  in  Fig.  28  summarizes  the  mean  adiabatic  MOC  circulation  for 
the  whole  ocean.  Fig.  29  is  a  schematic  of  the  inferred  zonally-averaged  meridional 
cross-section  of  the  whole  ocean. 
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Figure  27.  Adiabatic  interior  flux  Va  in  the  global  ocean  computed  from  the  four 
meteorological  databases. 
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The  Va  profiles  in  Fig.  27  shows  the  pattern  of  the  meridional  overturning  that  is 

generally  consistent,  in  terms  of  magnitude  and  direction,  with  earlier  numerical  and 
observational  studies. 

It  is  also  important  to  point  out  that  signs  and  magnitudes  of  Va  computed  from 
the  total  flux  and  their  wind-driven  components  based  on  the  Ekman  transport  are  poorly 
correlated.  One  of  the  interpretations  of  this  result  is  that  the  global  MOC  adiabatic 
component  may  not  be  mainly  driven  by  mechanical  forcing  generated  by  winds,  and  the 
eddy  contribution  is  essential  for  the  maintenance  of  meridional  overturning. 


Figure  28.  Schematic  of  the  adiabatic  meridional  overturning  circulation  in  the  global  ocean. 


The  Va  profile  computed  form  the  reference  climatology  (ECMWF)  indicates  that 

the  integrated  adiabatic  advection  over  the  whole  density  range  [cr(l,crmax  ]  is  almost  zero. 

This  implies  that,  as  expected,  the  northward  and  southward  interior  fluxes  balance  each 
other  in  this  density  interval,  lending  more  credence  to  the  proposed  diagnostics. 
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South 


North 


Figure  29.  Schematic  of  the  zonally-averaged  meridional  cross-section  of  the  global  ocean 
and  the  associated  MOC.  Red  arrows  label  southward  adiabatic  advection  and 
blue  arrows  depict  northward  adiabatic  advection. 


In  addition,  Fig.  27  highlights  the  discrepancies  among  the  Va  values  obtained 
from  the  four  climatologic  databases.  The  differences  in  the  sea  surface  density  flux 
between  these  datasets  do  not  result  in  as  large  changes  in  Va  profile  as  those  in  the  Vd 

computations.  Moreover,  we  still  observe  that,  similar  to  Vd  plots  (Fig.25),  the  Va  profiles 

based  on  ECMWF  and  NCEP/NCAR  climatology  are  relatively  close  to  each  other,  while 
those  calculated  from  SOC-GASC97  and  Da  Silva  SMD94  have  significant  differences. 

Table  4  gives  the  average  along-isopycnal  flux  differences  between  values 
obtained  from  ECMWF  re-analysis  and  those  computed  from  the  other  databases,  over 
the  potential  density  anomaly  interval  [24.2,  28.3].  As  in  section  IV-A-1,  results  in  this 
table  should  be  compared  with  those  in  Table  2.  This  comparison  shows  that  the 
estimates  of  isopycnal  volume  flux  are  much  less  dependent  on  the  density  flux  accuracy 
than  are  those  of  diapycnal  volume  flux.  Indeed,  the  Va  differences  among  the  different 

databases  are  much  lower  that  the  corresponding  Vd  differences.  Unlike  Vd ,  no  coherent 
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empirical  function  relating  the  absolute  mean  density  flux  accuracy  and  the  absolute 
average  adiabatic  flux  accuracy  was  found.  This  result  could  be  expected  since  signs  of 
mean  density  flux  differences  and  mean  isopycnal  volume  flux  differences  are  not 
correlated  [see  Tables  2  and  4], 


Mean  isopycnal  volume  flux  difference 

(mV) 

V A  ( NCEP /  NCAR)  -  V A  (ECMWF) 

-1.2506xl06  (*-l  Sv) 

VA  ( SOC  -  GASC97)  -  VA  (ECMWF) 

-4.0233x1 06  (  « -4  Sv) 

V A  (DaSilva)  -  V A  (ECMWF) 

-1.1483xl07  (*-11  Sv) 

Table  4.  Means  of  isopycnal  volume  flux  differences  between  ECMWF  re-analysis  and  the 
other  databases. 


3.  Adiabatic  Flux  Ratio 

Table  5  gives  the  adiabatic  flux  ratios  computed  from  each  meteorological 
climatology  using  Eq.  (28).  Here,  the  adiabatic  flux  ratios  from  ECMWF  and 
NCEP/NCAR  are  about  3.5  and  3.0,  respectively.  Their  magnitudes  support  the  theory 
that  the  global  Meriodional  Overturning  Circulation  has  a  dominant  adiabatic  component. 
The  SOC-GASC97  and  Da  Silva  SMD94  values  are  provided  here  only  for  references 
without  further  interpretation,  since  their  corresponding  Vd  and  Va  profdes  appear  to  be 
inconsistent  with  the  real  oceanic  circulation. 


Database 

Adiabatic  flux  ratio  y 
(dimensionless) 

ECMWF 

3.5018 

NCEP/NCAR 

3.0288 

SOC-GASC97 

1.4552 

Da  Silva  SMD94 

2.3247 

Table  5.  Adiabatic  flux  ratios  for  the  global  ocean  computed  from  the  four  meteorological 
databases. 
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B.  ATLANTIC  OCEAN 

If  we  limit  our  study  of  the  Atlantic  Ocean  to  the  density  range  for  which  all 
isopycnals  are  zonally  bounded  by  land  throughout  the  year,  the  potential  density 
anomaly  interval  reduces  to  [24.2,  25.4].  These  values  are  obtained  from  the  analysis  of 
the  monthly  mean  distribution  of  potential  density  [not  shown].  As  mentioned  previously, 
confining  the  analysis  to  this  density  range  permits  an  accurate  description  of  the  upper- 
ocean  meridional  circulation.  This  assures  the  absence  of  major  external  sources  of  water 
as  assumed  in  the  model  formulation  (section  II). 

Vd(cr)  and  Tj(cr)  over  the  sigma  range  [24.2,  28.3]  (used  for  the  global  ocean 

computations)  have  also  been  computed  in  order  to  assess  the  Atlantic  contribution  to  the 
global  MOC,  although  their  pattern  in  the  range  [26.5,  28]  could  be  affected  by  the 
contribution  from  the  Antarctic  Circumpolar  Current  (ACC). 

1.  Diapycnal  Volume  Flux 

The  diapycnal  interior  flux  ( Vd  )  in  the  upper  Atlantic  Ocean,  computed  from  the 

four  meteorological  databases  using  Eq.  (21)  and  Eq.  (22),  is  plotted  in  Fig.  30.  Its 
average  over  the  potential  density  anomaly  range  [24.2,  25.4]  is  approximately  6.7  Sv 
from  the  ECMWF  climatology  (reference  database). 

Fig.  30  shows  that  Vd  is  mostly  weak  over  the  density  interval  [24.20,  25.16]:  its 
values  remain  between  -9  Sv  and  +8  Sv  (ECMWF).  The  diabatic  volume  flux  abruptly 
increases  to  reach  its  maximum  magnitude  (37  Sv)  at  about cr  =  25.34  .  This  peak  value 
corresponds  to  the  high  upward  net  heat  flux  at  the  surface  of  the  Gulf  Stream.  Indeed, 
the  mean  potential  density  anomaly  at  40  m  depth  within  the  Gulf  Stream  is  cr «  25.35  . 
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Figure  30.  Diapycnal  interior  flux  Vd  in  the  Atlantic  Ocean  computed  from  the  four 
meteorological  databases. 


Similar  to  the  global  diapycnal  volume  flux  distribution,  we  note  large 
discrepancies  between  the  Vd  estimates  from  different  databases.  Results  from  ECMWF 
and  NCEP/NCAR  re-analyses  show  many  similarities,  but  NCEP/NCAR  density  forcing 
seems  to  be  too  high  with  regard  to  Vd  values.  Concerning  SOC-GASC97  and  Da  Silva 

SMD94  climatologies,  the  results  are  again  unrealistic,  which  we  attribute  to  the  lack  of 
observational  data  for  the  Southern  Atlantic  Ocean.  This  last  point  is  illustrated  in  Fig. 
31,  where  the  values  of  F(a)  =  jo  ■  /?  dl  [defined  in  Eq.  (13)]  in  the  Northern  and 

( 7 

Southern  Hemispheres  are  given  for  the  four  climatologic  databases:  the  differences  in 
Fn  values  are  much  lower  than  those  in  Fs  estimates. 
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2.  Adiabatic  Advection 

Fig.  32  gives  the  along-isopycnal  interior  fluxes  (Va),  computed  from  Eqs.  (21), 
(25),  and  (27)  using  the  four  meteorological  databases.  Here,  Va  (cr)  estimates  represent 
the  integrated  isopycnal  interior  flux  over  the  potential  density  anomaly  range  [cr0,cr] 
instead  of  over  the  density  interval  [cr,crmax]  (as  in  the  global  ocean  computations). 

Changes  in  the  direction  of  the  MOC  circulation  (northward  or  southward)  are 
still  observable  by  maxima  and  minima  in  the  Va  profiles,  but,  in  this  case,  the  direction 

dV  dV 

of  the  adiabatic  advection  at  each  cr  is  given  by  the  sign  of  — - :  — -  >  0  means  a 

da  da 

dV 

northward  flux,  and  -^-s-  <  0  corresponds  to  a  southward  flux.  Va  (cr) ,  calculated  from 

ECMWF  (reference  climatology),  shows  a  weak  southward  integrated  isopycnal  flow  for 
24.20  <  cr  <  25. 17  and  then  a  strong  northward  flow  for  25. 17  <  cr  <  25.38  .  This  points 
out  the  predominantly  northward  residual  flow  in  the  Atlantic  upper  thermocline.  Figure 
32  also  highlights  the  strong  effect  of  the  large  and  positive  density  flux  (into  the  ocean) 
in  the  Gulf  Stream  region  on  the  adiabatic  flow  in  the  upper  Atlantic  Ocean. 

Concerning  the  Ekman  fluxes,  the  plots  suggest  that  the  MOC  in  the  Atlantic 
Ocean  is  not  driven  by  the  local  mechanical  forcing  (wind  stress).  This  result  is  similar  to 
that  obtained  for  the  global  thennohaline  circulation. 
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Figure  31. 


Cross-isopycnal  volume  flux  within  the  mixed-layer,  (a)  In  the  Northern 
Hemisphere (Fn  (cr)) .  (b)  In  the  Southern  Hemisphere (  Fs  (ct))  . 
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Figure  32.  Adiabatic  interior  flux  Va  in  the  upper  Atlantic  Ocean  computed  from  the  four 
meteorological  databases. 


3.  Adiabatic  Flux  Ratio 

Table  6  gives  the  adiabatic  flux  ratios  for  the  Atlantic  basin.  This  parameter  is 
relatively  high  (y«3.6)  for  the  ECMWF  computations,  which  suggests  that  the 
meridional  overturning  in  the  Atlantic  thermocline  is  controlled  by  the  adiabatic 
processes  (advection  and  eddy-transfer).  However,  the  other  meteorological  databases  do 
not  support  this  hypothesis. 


Database 

Adiabatic  flux  ratio  y 
(dimensionless) 

ECMWF 

3.5675 

NCEP/NCAR 

1.1484 

SOC-GASC97 

0.3983 

Da  Silva  SMD94 

0.6259 

Table  6.  Adiabatic  flux  ratios  for  the  upper  Atlantic  Ocean  computed  from  the  four 
meteorological  databases. 
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4.  Atlantic  Contribution  to  Global  MOC 

In  order  to  estimate  the  Atlantic  contribution  to  the  global  thennohaline 
circulation,  we  can  compute  the  diabatic  and  adiabatic  volume  fluxes  over  the  density 
interval  [24.2,  28.3]  and  compare  them  to  the  values  corresponding  to  the  global  ocean. 

For  this  computation,  we  include  in  the  Atlantic  basin  the  southern  ocean  between 
Cape  Horn  and  Cape  of  Good  Hope.  Note  that  this  extension  of  the  Atlantic  surface  adds 
external  sources  of  water  and  might  lead  to  uncertainties  in  Vd  and  F  estimates.  In 

particular,  we  may  expect  differences  on  the  order  of  5-10  5V  between  the  entering  ACC 
flow  in  the  Drake  Passage  and  the  exiting  flow  to  the  south  of  Africa.  Such  differences 
are  not  significant  relative  to  the  ACC  mean  volume  flux  (100  5v),  but  may  be  relevant 
for  our  results  (on  the  order  of  10  5v). 

Figures  33  and  34  show  Vd  (a)  and  F  (a)  over  the  potential  density  anomaly 
range  of  [24.2,  28.3].  Note  that  Ffl(cr)  again  represents  the  MOC  adiabatic  advection 
over  the  density  range  [cr,crmax] ,  as  in  the  global  ocean  computations. 

Comparisons  between  Figs.  33  and  25  and  between  Figs.  34  and  27  show  many 
similarities  in  the  sign  and  magnitude  of  the  MOC  components.  These  similarities 
suggest  that  the  global  pole-to-pole  MOC  is  strongly  influenced  by  the  Atlantic  Ocean.  In 
particular,  the  peak  values  in  Vd  and  F  profiles  at  a  «  25.35 ,  which  correspond  to  the 
Gulf  Stream  signature,  are  noticeable  on  Figs.  25,  27,  33,  and  34. 

Table  7  provides  the  adiabatic  flux  ratios  for  the  entire  Atlantic  basin.  The  very 
high  value  (y  ~  4.11  for  ECMWF  climatology)  supports  the  conclusion  that  the 
meridional  overturning  in  the  whole  Atlantic  Ocean  is  mainly  adiabatic. 

Figures  32  and  34  display  the  discrepancies  among  the  Va  values  obtained  from 
the  four  climatologic  databases.  Again  we  see  that  the  differences  in  Va  (cr)  are  not  as 
large  as  those  observed  for  Vd  (cr) . 
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Figure  33.  Diapycnal  interior  flux  Vd  in  the  entire  Atlantic  Ocean  computed  from  the  four 
meteorological  databases. 
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Figure  34.  Adiabatic  interior  flux  Va  in  the  entire  Atlantic  Ocean  computed  from  the  four 
meteorological  databases. 
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Database 

Adiabatic  flux  ratio  y 
(dimensionless) 

ECMWF 

4.1055 

NCEP/NCAR 

2.5272 

SOC-GASC97 

1.7821 

Da  Silva  SMD94 

2.0053 

Table  7.  Adiabatic  flux  ratios  for  the  entire  Atlantic  Ocean  computed  from  the  four 
meteorological  databases. 


5.  Method  Comparison 

Section  III-A  presents  two  ways  for  computing  the  residual  flux  at  the  bottom  of 
the  mixed-layer.  The  first  method  uses  integration  along  the  isopycnal  contours  [Eq.  (14)] 
within  the  mixed-layer,  while  the  second  one  computes  integrals  over  the  surface 
bounded  by  these  contours  [Eq.  (21)]. 

The  two  methods  are  mathematically  equivalent,  but,  in  practice,  may  lead  to 
relatively  significant  differences  in  estimates  of  Va  and  Vd ,  depending  on  the  grid 

resolution.  Therefore,  comparing  their  results  gives  a  consistency  check  of  the  estimate 
accuracy.  Note  that  the  second  method  is  more  convenient,  and  for  this  reason  it  [Eq. 
(21)]  will  be  mostly  used. 

Figure  35  gives  Vd  (cr)  and  F  (a)  profiles  computed  by  the  two  methods  from 
ECMWF  climatology.  These  profiles  show  that  Va  and  Vd  estimates  weakly  depend  on 

the  method  used.  This  is  also  observed  for  the  adiabatic  flux  ratios,  which  are, 
respectively,  y  =  3.6607  and  y  =  3.5675,  for  the  integration  along  isopycnal  contours 
and  the  surface  integrals. 
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b) 


Figure  35.  MOC  components  in  the  Atlantic  Ocean  computed  from  ECMWF  climatology, 
(a)  Diapycnal  volume  flux  Vd  .  (b)  Adiabatic  volume  flux  Va . 


C.  PACIFIC  AND  INDIAN  OCEANS 

The  MOC  components  for  the  Indian  Ocean  have  not  been  computed  separately. 
Indeed,  the  Asian  continent  constitutes  a  northern  boundary  at  low  latitudes  and  therefore 
inhibits  any  isopycnal  transport  in  the  Indian  basin,  as  shown  by  the  distribution  of 
potential  density  at  40  m  depth  in  Fig.  9. 
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The  large  Pacific  Ocean  and  its  density  field  pennit  a  pole-to-pole  overturning 
circulation,  which  has  to  be  studied.  In  addition,  the  numerous  “holes”  in  its  western 
lateral  boundaries,  and  particularly  the  Indonesian  archipelago,  make  it  necessary  to  take 
into  account  the  density  forcing  over  the  Indian  Ocean  as  a  possible  source  of  water  mass 
transformation  in  the  Pacific  basin.  For  these  reasons,  calculations  will  be  carried  out  for 
the  Pacific  Ocean  only  and  for  the  combined  Indian-Pacific  Ocean,  and  results  compared. 

Similar  to  the  Atlantic  Ocean  computations,  we  first  limit  our  study  of  the  Pacific 
and  Indian  Oceans  to  the  density  range  for  which  all  isopycnals  are  zonally  bounded  by 
land  throughout  the  year:  the  monthly-averaged  potential  density  distribution  [not  shown] 
provides  the  potential  density  anomaly  range  of  [23.8,  25.4].  It  is  also  interesting  to 
assess  the  Pacific  and  Indian  contributions  to  the  global  MOC  by  computing  Vd  (cr)  and 

Va  (cr)  over  the  sigma  range  [24.2,  28.3],  used  for  the  global  ocean  calculations. 

1.  Diapycnal  Volume  Flux 

The  diapycnal  interior  flux  (  Vd  )  in  the  upper  Pacific  Ocean  is  plotted  in  Fig.  36. 

Its  average  value  over  the  potential  density  anomaly  range  [23.8,  25.4]  is  about  3.6  Sv 
(ECMWF  climatology).  This  mean  value  is  less  than  that  computed  in  the  upper  Atlantic 
Ocean,  but,  at  each  cr ,  Vd  is  generally  higher  in  the  Pacific  than  in  the  Atlantic. 

We  note  some  peak  values  on  the  Vd  profile,  but  no  characteristic  signature  of  the 

Kuroshio  Current  (cr«  25.0).  Indeed,  we  would  have  expected  a  positively  high  value 

around  a  =  25.0  due  to  the  high  upward  net  heat  flux  at  the  surface  of  this  warm  western 
boundary  current  [see  Fig.  10].  Instead,  we  observe  a  negative  peak  value  at  cr «  25.1. 
This  result  suggests  that  the  Gulf  Stream  has  a  much  bigger  impact  on  the  MOC  than  the 
Kuroshio  Current.  In  fact,  Fig.  10  indicates  that  the  highest  density  forcing  is  in  the  Gulf 
Stream  region. 
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Figure  36.  Diapycnal  interior  flux  Vd  in  the  Pacific  Ocean  computed  from  the  four 
meteorological  databases. 


Again,  we  observe  large  discrepancies  among  ^estimates  from  the  different 

databases,  although  results  from  ECMWF  and  NCEP/NCAR  re-analyses  show  many 
similarities.  Concerning  SOC-GASC97  and  Da  Silva  SMD94  climatologies,  the  results 
are  again  unrealistic,  which  is  most  likely  due  to  the  lack  of  observational  data  for  the 
southern  ocean. 

The  Vd  profiles  for  the  combined  Indian-Pacific  Ocean  are  close  to  those  for  the 

Pacific  Ocean  only.  As  shown  in  Fig.  37,  the  major  difference  is  in  the  value  of  the 
diapycnal  interior  flux,  which  increases  by  13  Sv  (ECMWF  climatology)  when  the  Indian 
Ocean  is  taken  into  account.  This  implies  that  the  Indian  Ocean  contribution  to  Vd  (cr)  is 
almost  constant  (about  13  Sv)  over  the  sigma  range  [23.8,  25.4]. 
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Figure  37. 


Comparison  of  the  diapycnal  interior  flux  Vd  in  the  Pacific  Ocean  and  in  the 


combined  Indian-Pacific  Ocean.  The  values  are  computed  from  ECMWF  and 
NCEP/NCAR  databases. 


2.  Adiabatic  Advection 

Figure  38  shows  Fa(cr)  computed  from  Eqs.  (21),  (25),  and  (27)  using  the  four 
meteorological  databases.  In  a  similar  way  as  in  the  upper  Atlantic  calculations,  Va  (cr) 

represents  here  the  integrated  isopycnal  interior  flux  over  the  potential  density  anomaly 

,  ,  dV  dV 

range  I  cr0,cr  I .  Thus,  — -  >  0  means  a  northward  flux,  and  — -  <  0  corresponds  to  a 
da  da 

southward  flux. 

The  integrated  isopycnal  flux  Va,  calculated  from  ECMWF  (reference 
climatology),  is  negative  over  most  of  the  cr  range  [23.8,  25.4].  This  indicates  that  the 
residual  flow  in  the  Pacific  upper  thennocline  is  predominantly  southward,  with  a  mean 
value  of  -8.6  Sv.  The  differences  in  Va  (cr)  computed  from  the  various  databases  are  not 

as  large  as  in  Vd  (cr) ,  which  is  consistent  with  our  results  for  the  Atlantic  Ocean. 
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Contrary  to  the  results  in  the  Atlantic  and  global  oceans,  the  comparison  of  the 
ispoycnal  Ekman  fluxes  and  the  isopycnal  residual  (total)  fluxes  suggests  that  mechanical 
forcing  (wind  stress)  is  the  major  driving  force  for  the  MOC  in  the  Pacific  Ocean.  Indeed, 
the  Ekman  flux  dominates  the  residual  flux,  with  the  same  sign  and  tendency,  over  almost 
the  complete  density  range.  Note  that  this  statement  is  not  valid  for  the  combined  Indian- 
Pacific  basin. 

The  computations  for  the  combined  Indian-Pacific  Ocean  provide  Va  profiles  that 
are  similar  to  those  for  the  Pacific  Ocean  only.  As  shown  in  Fig.  39,  we  observe  a  nearly 
uniform  Indian  Ocean  contribution  to  Va  (cr)  that  is  about  12  Sv  over  the  entire  sigma 
range  [23.8,  25.4]. 
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Figure  38.  Adiabatic  interior  flux  Va  in  the  upper  Pacific  Ocean  computed  from  the  four 
meteorological  databases. 
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Figure  39.  Comparison  of  the  adiabatic  interior  flux  Va  in  the  Pacific  Ocean  and  in  the 

combined  Indian-Pacific  Ocean.  The  values  are  computed  from  ECMWF  and 
NCEP/NCAR  databases. 


3.  Adiabatic  Flux  Ratio 

Tables  8  and  9  give  the  adiabatic  flux  ratios  respectively  for  the  Pacific  basin  and 
the  combined  Indian-Pacific  Ocean.  This  parameter  is  relatively  low  in  the  two  oceans: 
y~2A5  (ECMWF  computations)  in  only  the  Pacific  Ocean  and  /«1.47  (ECMWF 
computations)  when  we  also  consider  the  Indian  basin.  This  suggests  that  the  meridional 
overturning  is  still  controlled  by  the  adiabatic  processes  (advection  and  eddy-transfer)  in 
the  Pacific  thermocline,  whereas  it  is  mostly  diabatic  in  the  upper  Indian  Ocean. 
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Database 

Adiabatic  flux  ratio  y 
(dimensionless) 

ECMWF 

2.4542 

NCEP/NCAR 

1.4346 

SOC-GASC97 

0.4461 

Da  Silva  SMD94 

0.5854 

Table  8.  Adiabatic  flux  ratios  for  the  upper  Pacific  Ocean  computed  from  the  four 
meteorological  databases. 


Database 

Adiabatic  flux  ratio  y 
(dimensionless) 

ECMWF 

1.4681 

NCEP/NCAR 

0.7260 

SOC-GASC97 

0.5042 

Da  Silva  SMD94 

0.5683 

Table  9.  Adiabatic  flux  ratios  for  the  upper  Indian-Pacific  Ocean  computed  from  the  four 
meteorological  databases. 


4.  Pacific  and  Indian  Contributions  to  Global  MOC 

To  estimate  the  Pacific  and  Indian  contributions  to  the  global  thennohaline 
circulation,  Figs.  40,  41,  42,  and  43  display  Vd  (cr)  and  Va  (cr)  in  the  Pacific  and  Indian- 

Pacific  oceans  over  the  density  interval  of  [24.2,  28.3].  The  chosen  density  range  includes 
the  southern  ocean,  where  the  volume  flux  of  the  ACC  may  result  in  external  sources  of 
water  and  might  lead  to  uncertainties  in  Vd  and  Va  estimates. 
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Figure  40. 


Diapycnal  interior  flux  Vd  in  the  entire  Pacific  Ocean  computed  from  the  four 
meteorological  databases. 
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Figure  41.  Adiabatic  interior  flux  Va  in  the  entire  Pacific  Ocean  computed  from  the  four 
meteorological  databases. 
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In  Figs.  41  and  43,  Va(cr)  represents  the  integrated  adiabatic  advection  over  the 
density  range  [cr,crmax].  In  Fig.  41,  we  see  that  Va  is  dominated  by  the  Ekman  flux. 
Significant  differences  between  Va  and  Va_EK  occur  only  for  the  low  densities  (cr  <  25.3) : 
this  feature  may  be  partially  due  to  the  external  water  sources  in  the  southern  ocean 
affecting  the  integrated  flow  Va  (cr) .  The  mechanical  forcing  (wind  stress)  also  drives  the 

MOC  in  the  combined  Indian-Pacific  Ocean,  except  at  low  densities  (cr  <  25.2) .  The 

Ekman  flux  dominates  the  residual  flux,  with  the  same  tendency,  over  almost  all  the 
density  range  [25.2,  28.3], 

The  comparisons  between  Fig.  42  and  Fig.  25  and  between  Fig.  43  and  Fig.  27 
indicate  the  large  contribution  of  the  density  forcing  over  the  Indian  Ocean  to  the  global 
MOC  circulation  (on  both  diabatic  and  adiabatic  components).  Indeed,  the  correlations 
between  the  global  MOC  components  and  those  computed  over  the  Pacific  Ocean  only 
are  relatively  low,  whereas  the  same  correlations  with  the  combined  Indian-Pacific  Ocean 
instead  of  only  the  Pacific  Ocean  are  much  higher.  This  result  is  surprising,  since  the 
Indian  Ocean  is  the  smallest  oceanic  basin  and  its  topography  also  inhibits  any  pole-to- 
pole  MOC  circulation. 

Tables  8  and  9  give  the  adiabatic  flux  ratios  respectively  for  the  whole  Pacific 
basin  and  the  whole  combined  Indian-Pacific  Ocean.  Relatively  high  values  of  y  indicate 
that  the  adiabatic  advection  is  also  an  important  process  in  the  Pacific  and  Indian  Oceans. 
In  particular,  the  isopycnal  volume  flux  tends  to  dominate  the  MOC  at  high  densities 
(southern  ocean). 
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Figure  42.  Comparison  of  the  diapycnal  interior  flux  Vd  in  the  entire  Pacific  Ocean  and  in 

the  combined  Indian-Pacific  Ocean.  The  values  are  computed  from  ECMWF  and 
NCEP/NCAR  databases. 


Figure  43.  Same  as  Fig.  42,  but  for  the  adiabatic  interior  flux  Va . 
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Database 

Adiabatic  flux  ratio  y 
(dimensionless) 

ECMWF 

2.5561 

NCEP/NCAR 

2.8388 

SOC-GASC97 

0.9889 

Da  Silva  SMD94 

1.6378 

Table  10.  Adiabatic  flux  ratios  for  the  entire  Pacific  Ocean  computed  from  the  four 
meteorological  databases. 


Database 

Adiabatic  flux  ratio  y 
(dimensionless) 

ECMWF 

3.0824 

NCEP/NCAR 

3.1211 

SOC-GASC97 

1.5104 

Da  Silva  SMD94 

2.3517 

Table  11.  Adiabatic  flux  ratios  for  the  entire  Indian-Pacific  Ocean  computed  from  the  four 
meteorological  databases. 


D.  CONSISTENCY  OF  THE  ADIABATIC  ADVECTION  ESTIMATES 

The  isopycnal  flux  Va  (cr  )  is  defined  by  a  linear  average  of  two  net  residual  fluxes 

at  the  bottom  of  the  mixed-layer  [Eq.  (25)].  Since  this  assumption  is  a  critical  element  of 
our  method  for  estimating  the  residual  overturning  circulation,  the  formulation  for  Va 
must  be  tested  for  consistency. 

One  way  to  test  Va  estimates  is  to  use  an  oceanic  General  Circulation  Model 
(GCM).  Indeed,  the  residual  isopycnal  advection  diagnosed  from  the  air-sea  density 
fluxes  can  be  compared  with  the  Va  values  derived  directly  from  the  GCM-simulated 
velocities. 


1.  Model  Characteristics 

The  model  used  for  this  study  is  based  on  the  GFDL  MOM3  code  (Pacanowski 


and  Griffies  1999).  Its  horizontal  resolution  is  two  degrees  in  longitude  and  latitude.  It 
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covers  the  whole  ocean  from  78°S  to  84°N  and  is  divided  into  25  levels  in  the  vertical 
direction,  with  its  resolution  increasing  from  17  m  at  the  surface  layer  to  510  m  at  the 
bottom.  The  bathymetry  of  the  model  is  derived  from  the  Scripps  Topography.  Vertical 
diffusivity  varies  from  0.25x10  4  nr.s  1  at  the  surface  to  1 .0x10  4  m2s  1  at  the  bottom. 
This  profile  reflects  the  increase  of  the  vertical  mixing  from  the  thermocline  to  the  deep 
ocean  (Bryan  and  Lewis  1979)  and  the  intensification  of  mixing  by  rough  bottom 
topography  (Polzin  et  al.  1997).  Heat  and  salt  transports  by  mesoscale  eddies  are 
parameterized  by  the  Gent-McWilliams  scheme  (Gent  and  McWilliams  1990)  with 
coefficients  for  isopycnal  diffusion  of  tracers  and  isopycnal  thickness  of  500  ms _1 .  The 
K-profile  parameterization  (KPP)  scheme  (Large  et  al.  1994)  is  used  to  represent 
turbulent  mixing  within  boundary  layers. 

The  simulated  density  field  at  40  m  depth  is  given  in  Fig.  44a:  due  to  the  coarse 
resolution  of  the  model,  small  differences  with  observations  remain  at  high  latitudes  of 
the  North  Atlantic  and  in  the  ice-covered  areas.  Surface  density  fluxes  are  quantitatively 
consistent  with  observations,  and  the  largest  differences  with  the  observational  data 
appear  in  regions  where  the  coarse  resolution  distorts  simulations  of  the  currents  (at  the 
equator,  in  the  Western  Boundary  Current  areas,  and  in  the  Antarctic  Circumpolar 
Current  region). 

Heat  fluxes  into  the  ocean  are  determined  by  using  conventional  bulk  formulas 
and  are  corrected  for  stability.  Daily  values  for  the  two-meter  air  temperature  and 
humidity,  ten-meter  wind  speed,  and  zonal  and  meridional  components  of  the  wind  stress 
have  been  taken  from  NCEP/NCAR  reanalysis  for  the  period  of  1979-2001  and 
climatologic  monthly  means  are  used  for  all  other  atmospheric  variables.  Cloud  cover 
and  solar  radiation  are  given  by  the  International  Satellite  Cloud  Climatology  Project 
(ISCCP).  Freshwater  fluxes  are  taken  from  Jiang  et  al.  (1999)  and  include  river  runoff 
data.  The  model  is  coupled  to  a  thermodynamic  sea-ice  model  (Visbeck  et  al.  1998). 

For  the  analysis,  the  1990  simulation  has  been  chosen:  it  represents  the  mean  year 
of  the  studied  period  and  does  not  correspond  to  an  El  Nino/La  Nina  event.  The  density 
flux  distribution  computed  from  the  model  is  given  in  Fig.  44b. 
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The  MOC  estimates  provided  by  the  model  are  generally  coherent  with  the 
previous  studies.  In  addition,  density  fluxes,  stratification,  and  velocity  fields  in  GCM- 
simulation  are  all  dynamically  consistent  and  readily  available,  which  enables  us  to  verify 
main  assumptions  of  the  theory  in  Section  II. 

2.  Adiabatic  Advection 

The  isopycnal  flux  has  been  estimated  over  the  interval  [24.4,  28],  This  range  is 
determined  with  the  same  criteria  as  those  used  in  Section  III-A-3:  a0  =24.4  is  the 
lightest  density  anomaly  whose  outcrop  forms  two  non-intersecting  contours  of  surface 
density  throughout  the  year,  and  crmax  =28  is  the  highest  density  of  the  ice-free  sea 
surface  area. 

The  adiabatic  MOC  component  Va  (cr)  has  been  computed  in  two  ways:  (i)  by 
using  Eq.  (26)  to  infer  Va  from  the  sea-surface  density  flux,  and  (ii)  by  evaluating  Va 

directly  from  the  meridional  velocity  distribution.  More  precisely,  for  the  second 
calculation,  we  consider  a  zonal  section  exactly  at  the  equator,  and  the  isopycnal 
transport  is  obtained  by  vertically  integrating  the  zonal  residual  flow  (the  sum  of  the 
Eulerian  and  eddy-induced  velocities  of  the  Gent-McWilliams  scheme)  starting  from 
<7max  towards  lighter  densities.  The  sign  convention  in  both  calculations  is  the  same  as  in 
Section  IV-A-2. 
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Figure  44.  GCM-simulated  fields,  (a)  Potential  density  anomaly  distribution  at  40  m  depth. 

(b)  Surface  density  flux  from  the  atmosphere  into  the  ocean  (from  Radko  et  al. 
2006). 
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For  the  individual  basin  computations,  the  water  mass  transport  from  other  basins 
has  been  taken  into  account.  In  particular,  the  ACC  contribution  has  been  calculated  as  a 
divergence  of  the  zonal  residual  flow  within  the  ACC  between  64°W  and  18°E,  which  are 
the  longitudes  of,  respectively,  the  Drake  Passage  and  Africa.  Over  the  southern  Atlantic 
area,  we  observe  a  negative  divergence:  this  suggests  some  inflow  of  water  from  the 
ACC  into  the  Atlantic  Ocean.  This  ACC  convergence  has  been  incorporated  into  the  Vs 
tenn  of  Eq.  (25)  for  the  Atlantic  and  Indian-Pacific  calculations,  and  the  corresponding 
Va  estimate  has  been  modified  accordingly. 

Figures  45,  46,  and  47  plot  Va  (cr)  for  the  global  ocean,  the  Atlantic  and  Indian- 
Pacific  basins,  respectively.  On  these  plots,  we  see  that  the  Va  profiles  from  the  surface 

density  flux  (blue  curves)  and  those  obtained  directly  from  the  velocity  distribution  at  the 
equator  (green  curves)  match  very  well,  especially  for  the  global  calculations.  Although 
some  significant  differences  occur  at  specific  cr  values,  overall,  the  results  in  Figs.  45, 
46,  and  47  support  the  diagnostic  model  for  estimating  Va  (cr) ,  introduced  in  Section  II-D 
[Eq.  (25)]. 

For  the  global  ocean  case,  the  profile  of  the  adiabatic  MOC  component  computed 
from  observations  [Fig.  27]  and  that  based  on  model  calculations  [Fig.  45]  have  many 
similarities.  In  particular,  we  observe  on  both  profiles  a  positive  peak  value  of  5  Sv  at 
cr  «  25.7  and  a  negative  peak  value  at  cr  ~  27. 1  in  Fig.  27  (observations)  and  at  cr  ~  27.3 
in  Fig.  45  (model).  However,  the  magnitude  of  Va  (cr)  computed  from  observations  is 
generally  higher  than  in  the  model.  The  Va  profile  given  by  the  model  does  not  show  any 
Gulf  Stream  signature  (negative  peak  value  at  cr  «  25.35  ),  whereas  it  clearly  appears  on 
the  observational  Va  curve.  This  may  be  attributed  to  poor  resolution  of  WBC  Gulf 
Stream  density  field  by  the  coarse-resolution  GCM. 

For  the  Atlantic  Ocean  computations,  Efl  (cr)  profiles  from  observations  [Fig.  35] 

and  GCM  [Fig.  46]  are  quite  similar,  but,  as  mentioned  above,  they  strongly  differ  at  the 
Gulf  Stream  densities  (around  cr  =  25.35  ). 
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Va  (cr)  profiles  for  the  Indian-Pacific  Ocean  [Fig.  43  and  Fig.  47]  also  show  many 

similarities  between  the  observations  and  the  model  output.  Again,  the  magnitude  is 
generally  much  higher  for  the  climatologic  estimates  than  for  the  GCM-based  analysis. 

Concerning  the  wind-driven  component  of  Va  (cr)  (Ekman  isopycnal  flux),  the 

model  outputs  are  characterized  by  a  much  larger  contribution  of  this  component  to  the 
residual  (total)  isopycnal  flux  than  the  observational  estimates.  Flowever,  both  the 
observational  results  and  the  model  suggest  that  the  local  wind-stress  forcing  does  not 
drive  the  residual  MOC  adiabatic  advection  in  the  Atlantic  Ocean. 


Global 


Figure  45.  Adiabatic  interior  flux  Va  in  the  global  ocean  computed  from  the  GCM.  The 
displayed  values  are  in  sverdrups  (SV).  The  blue  curve  represents  Va  (cr) 
computed  from  the  simulated  surface  density  flux.  The  green  curve  gives 
Va  (cr)  obtained  directly  from  the  velocity  distribution  at  the  equator.  The  red 

curve  is  the  wind-driven  component  of  Va  (cr)  (Ekman  adiabatic  interior  flux) 
[from  Radko  et  al.  2006]. 
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Atlantic 


Figure  46.  Same  as  Fig.  45,  but  for  the  Atlantic  Ocean. 


Indo-Pacific 


Figure  47.  Same  as  Fig.  45,  but  for  the  Indian-Pacific  Ocean. 
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V.  DISCUSSION  AND  CONCLUSIONS 


The  oceanic  Meridional  Overturning  Circulation  (MOC)  is  accompanied  by  the 
diabatic  transformation  of  water  masses  and  interacts  with  the  air-sea  buoyancy  fluxes 
(Walin  1982).  The  two  major  mechanisms  for  water  mass  transformation  are  the  interior 
diapycnal  diffusion  and  the  diabatic  processes  in  the  upper  mixed-layer.  Their  relative 
importance  is  still  uncertain  and  much  debated.  The  classical  view  (e.g.,  Robinson  and 
Stommel  1959;  Munk  1966)  emphasizes  the  water  mass  transfonnation  in  the  ocean 
interior  and  implicitly  assigns  it  a  critical  role  in  maintenance  of  the  meridional 
overturning.  This  idea  was  challenged  by  a  number  of  modeling  studies  (e.g.,  Bryan 
1987),  which  demonstrate  that  reproducing  overturning  of  realistic  strength  via  diffusive 
mechanisms  requires  unrealistically  high  values  of  vertical  mixing.  Mechanical  forcing 
by  winds,  on  the  other  hand,  appears  to  be  an  essential  ingredient  of  the  thermohaline 
system;  neglecting  the  wind  stress  in  multi-century  simulations,  for  instance,  leads  to  a 
complete  shutdown  of  the  conveyor  belt  circulation  (Timmermann  and  Goosse  2004). 
These  findings  suggest  that  the  interior  diffusion  may  not  be  as  critical  for  maintenance 
of  the  MOC  as  adiabatic  processes  (wind  forcing,  advection,  and  eddy  transfer).  This 
latter  statement  has  gained  considerable  attention  recently  (Toggvieler  and  Samuels  1998; 
Marshall  and  Radko  2003;  Radko  2006),  but  requires  independent  and  objective  support. 

The  diagnostic  analysis  in  this  study  may  help  to  resolve  this  dispute  by 
quantifying  the  magnitude  and  distribution  of  the  isopycnal  MOC  component  (Va)  and 

relating  it  to  the  surface  forcing  fields.  While  the  isopycnal  transport  is  not  reflected  in 
the  interior  water  mass  transformation,  it  interacts  with  the  diabatic  mixed-layer,  making 
it  possible  to  infer  its  pattern  directly  from  the  sea  surface  data.  The  theoretical 
framework  is  based  on  the  residual  mean  theory  that  takes  into  account  adiabatic  transfer 
by  the  geostrophic  eddies.  As  indicated  in  the  schematic  diagram  in  Fig.  6b,  each  density 
layer  in  our  model  can  be  thought  of  as  a  “leaky  pipe”  that  is  connecting  the  two 
hemispheres  and  is  pumped  at  its  two  ends  with  different  forces,  ft  is  this  difference  in 
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pumping  that  ultimately  determines  the  direction  and  the  average  transport  of  the  flow 
through  the  pipe,  although  a  certain  amount  of  water  “leaks”  in  the  interior  due  to 
diabatic  processes. 

This  study  is  mainly  focused  on  the  along-isopycnal  volume  flux  Va ,  representing 
the  adiabatic  advection  within  a  specific  density  layer.  Va  (cr)  has  been  defined  as  the 

linear  average  of  the  two  net  residual  fluxes  at  the  bottom  of  the  mixed-layer  in  the 
Northern  and  Southern  hemispheres.  Moreover,  considerable  attention  is  given  to 
properties  of  the  diapycnal  flux  Vd(cr),  representing  the  integrated  flux  across  the 

isopycnal  surface  a  in  the  ocean  interior.  Since  the  residual  flow  is  assumed  to  be  non- 
divergent,  it  is  exactly  equal  to  the  net  volume  flux  from  the  ocean  interior  into  the 
mixed-layer  within  the  area  defined  by  the  mixed-layer  density  condition  am  <  <j  . 

In  addition  to  the  computations  of  the  adiabatic  and  diabatic  interior  fluxes,  the 
adiabatic  flux  ratio  defined  in  Eq.  (28)  has  been  used  to  determine  the  dominant 
component  of  the  MOC.  This  parameter  appears  to  represent  accurately  the  relative 
strength  of  the  two  components  by  considering  the  averages  of  the  volume  fluxes. 
However,  it  is  very  sensitive  to  the  accuracy  of  the  observational  fields  and  also  to  the 
considered  density  range.  For  this  reason,  the  values  of  the  adiabatic  flux  ratio  y  must  be 
interpreted  with  caution. 

For  the  Atlantic,  our  results  show  that  along-isopycnal  volume  fluxes  greatly 
exceed  the  corresponding  diapycnal  transport  (y  «  4.1 1  for  ECMWF  climatology).  These 
findings  support  the  notion  that  the  dynamics  of  the  meridional  overturning  in  the 
Atlantic  Ocean  are  dominated  by  the  adiabatic  processes  (advection  and  eddy-transfer). 
For  the  Pacific,  the  analysis  suggests  that  the  effects  of  the  water  mass  transformation  and 
along-isopycnal  advection  are  mostly  comparable  (y  ~  2.56  for  ECMWF  climatology). 
For  the  global  ocean,  the  results  suggest  that  the  global  thermohaline  circulation  is 
dominated  by  the  adiabatic  advection  (y  ~  3.50  for  ECMWF  climatology). 

In  addition  to  quantifying  the  strengths  of  the  MOC  components,  this  study  offers 
some  insight  into  the  mechanics  of  the  pole-to-pole  meridional  overturning.  The  total 
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along-isopycnal  residual  circulation  (Va)  is  divided  into  two  distinct  components:  the 
mean  advection  (Va_EK),  also  called  the  wind-driven  component  of  Va,  and  the  eddy- 
driven  component.  The  GCM  results  point  out  that  Va  is  controlled  by  the  mean 

advection,  directly  associated  with  the  Ekman  pumping,  everywhere  except  in  the 
Atlantic  Ocean.  On  the  contrary,  the  observational  results  show  that  the  Ekman  advection 
is  comparatively  weak  and  not  correlated  with  the  residual  adiabatic  advection,  except  in 
the  Pacific  Ocean;  this  emphasizes  the  global  dominance  of  the  eddy-driven  advection  in 
the  adiabatic  interior  flux. 

It  is  also  essential  to  keep  in  mind  that  the  application  of  the  water  mass 
transfonnation  theory  to  the  sea  surface  data  is  greatly  restricted  by  the  observational 
uncertainties:  the  air-sea  fluxes  in  various  datasets  may  differ  by  as  much  as  a  factor  of 
two  or  more.  However,  this  limitation  is  not  nearly  as  severe  for  the  isopycnal  MOC 
component  ( Va  )  as  it  is  for  the  diapycnal  flux  ( Vd ).  Va  estimates  are  based  on  the 
difference  between  the  air-sea  fluxes  at  the  Northern  and  Southern  outcrops  of  each 
isopycnal  surface.  Therefore,  for  Va(cr)  computations,  our  method  results  in  a  partial 

cancellation  of  systematic  errors  in  flux  measurements.  In  fact,  the  four  climatologic 
databases  (ECMWF,  NCEP/NCAR,  SOC-GASC97,  and  Da  Silva  SMD94)  are  mutually 
consistent  in  tenns  of  predicting  the  isopycnal  transport  Va . 

Several  databases  have  been  used  in  this  study:  Levitus  94  for  the  oceanic 
parameters,  ECMWF,  NCEP/NCAR,  SOC-GASC97,  and  Da  Silva  SMD94  for  the  air-sea 
fluxes.  This  has  pennitted  a  quality  control  of  the  results  and  an  assessment  of  the 
sensitivity  of  the  calculations  to  the  observational  fields.  Furthermore,  the  computations 
based  on  GCM-simulated  parameters  has  validated  our  method  of  estimating  the 
adiabatic  interior  flux  (Va  (cr)j . 

Statistical  calculations  in  this  study  offer  suggestions  for  the  database  selection 
and  for  the  optimal  choice  of  the  operational  parameters  in  diagnostic  models.  In 
particular,  the  mixed-layer  density  here  is  represented  by  the  potential  density  anomaly  at 
40  m  depth,  and  ECMWF  re-analysis  appears  to  be  the  best  compromise  regarding  the 
accuracy  of  fluxes  at  the  sea  surface. 
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Finally,  this  study  offers  new  insights  into  mechanics  of  the  thennohaline 
circulation  and  its  interaction  with  the  atmosphere.  These  insights  are  supported  by  an 
agreement  (Section  IV)  between  our  results  and  the  corresponding  estimates  of  the 
oceanic  circulation  from  in-situ  measurements  and  from  the  earlier  Walin-type  diagnostic 
calculations  (Donners  et  al.  2005). 
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